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The symmetry of magnetic quantum tunneling (MQT) in the single molecule magnet Mni2-acetate 
has been determined by sensitive low-temperature magnetic measurements in the pure quantum 
tunneling regime and high frequency EPR spectroscopy in the presence of large transverse magnetic 
fields. The combined data set definitely establishes the transverse anisotropy terms responsible 
for the low temperature quantum dynamics. MQT is due to a disorder induced locally varying 
quadratic transverse anisotropy associated with rhombic distortions in the molecular environment 
(2"'* order in the spin-operators). This is superimposed on a 4 order transverse magnetic anisotropy 
consistent with the global (average) S4 molecule site symmetry. These forms of the transverse 
anisotropy are incommensurate, leading to a complex interplay between local and global symmetries, 
the consequences of which are analyzed in detail. The resulting model explains: (1) the observation 
of a twofold symmetry of MQT as a function of the angle of the transverse magnetic field when a 
subset of molecules in a single crystal are studied; (2) the non-monotonic dependence of the tunneling 
probability on the magnitude of the transverse magnetic field, which is ascribed to an interference 
(Berry phase) effect; and (3) the angular dependence of EPR absorption peaks, including the fine 
structure in the peaks, among many other phenomena. This work also establishes the magnitude 
of the 2"** and 4*'' order transverse anisotropy terms for Mni2-acetate single crystals and the angle 
between the hard magnetic anisotropy axes of these terms. EPR as a function of the angle of the field 
with respect to the easy axis (close to the hard-medium plane) confirms that there are discrete tilts 
of the molecular magnetic easy axis from the global (average) easy axis of a crystal, also associated 
with solvent disorder. The latter observation provides a very plausible explanation for the lack of 
MQT selection rules, which has been a puzzle for many years. 

PACS numbers: 75.45.-|-j, 75.50.Tt, 75.60.Lr 
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I. INTRODUCTION 

The origin of magnetic quantum tunneling (MQT) in 
single-molecule magnets (SMMs) is important for funda- 
mental reasons as well as proposed applications, ranging 
from magnetic data storage to quantum computing 
A macroscopic (millimeter-sized) SMM single crystal to 
a first approximation can be considered an ensemble 
of weakly interacting molecules with the same chemical 
composition and orientation. In this regard SMM single 
crystals are ideal materials in which to study the quan- 
tum properties of magnetic molecules, as the magnetic 
response is amplified by the huge number of molecules 
forming the crystal. The environment of the molecules 
in crystals is also well defined and can be characterized 
by techniques such as x-ray diffraction. 

SMMs have a predominate uniaxial magnetic 
anisotropy that determines an easy magnetic axis 



for the spin. MQT is due to interactions that break 
this axial symmetry and lead to transitions between 
magnetic states with opposite projections of spin on 
the magnetic easy axis. While this is fundamental 
to their quantum dynamics, remarkably only recently 
have the nature of the transverse interactions that 
produce MQT been determined in the first and most 



Small modulations in the local environment around the 
magnetic cores have been found to be important in MQT 
(16. 17, 1^ llSl]. In particular, it has also been shown 
that in Mni2-ac there are a variety of types of disorder 
that affect the magnetic properties of the molecules, 
including a distribution of solvent microenvironments, 
g~ and D-strain ^HllillllllllilHI and a 
distribution of tilts of the easy axes of the molecules 
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FIG. 1: Energy diagram of the 25* + 1 = 21 projections of 
the spin 5 = 10 along the easy magnetic axis of the Mni2- 
ac molecule for resonances fc = and k = 1. The tunnel 
sphtting, A„ „j/ (m + m' = fc), is illustrated in both cases 
(not to scale). The inset shows the arrangement of Mn ions 
looking down the S4 axis of the molecule (z-axis). Note that 
the spins point up and down along this axis, perpendicular to 
the plane of the page. 



In this paper we present high sensitivity magnetome- 
try and high frequency electron paramagnetic resonance 
(EPR) studies carried out on single crystals of Mni2-ac. 
These experimental techniques and an in-depth analysis 
of the implications of the combined data set allow us to 
determine the symmetry and origin of MQT. We show 
that variations in the local solvent environments around 
the Mni20i2 magnetic cores are important to MQT be- 
cause they lower the symmetry of the core and lead to 
tilts of the magnetic easy axis. Disorder in the molecule's 
solvent environment can explain many of the features 
observed experimentally. This includes non exponential 
magnetic relaxation, the absence of tunneling selection 
rules, an unusual Berry phase effect, and multiple and 
broad EPR absorption peaks. 

The article is organized as follows: In Section II we 
analyze the spin Hamiltonian of Mni2-ac and discuss the 
symmetry of MQT expected based on this Hamiltonian. 
The effect of a disorder induced transverse anisotropy 
combined with intrinsic transverse anisotropy has inter- 
esting consequences for MQT that we explore in depth. 
In Section III we describe magnetic relaxation measure- 
ments carried out in the pure quantum regirne, where the 
MQT is not assisted by thermal activation [ll|. Section 
IV presents high frequency EPR measurements. In Sec- 
tion V we discuss the implications of our analysis and 
measurements to the understanding of MQT. 



II. 



TUNNELING IN MN12-AC 



Mni2-ac consist of a core of four Mn^"* ions {S = 3/2) 
surrounded by a ring of eight Mn^'^ ions (spin S = 2) (see 
fig. 1). These two groups of spins order ferrimagnetically 
producing: a net spin S = 10 (8x2-4x3/2 = 10) 

The spin Hamiltonian of [Mni20i2(CH3COO)i6 
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FIG. 2: Magnetization versus field for a single crystal of Mni2- 
ac for several sweep rates, a — dH/dt. The measurements 
have been conducted in the pure quantum regime at T = 
0.6 Kin. The observed increases of the magnetization at 
regularly spaced field intervals, Hk — D/g^B ~ 0.44 T, cor- 
respond to MQT at resonances fc = 6, 7, 8, 9 and 10. 



(H20)4]-2CH3COOH-4H20 is given by 

n = -DSl - BSt - g^JLB^ioH,S, cos d + HT + HA + n' , 

(1) 

The first two terms represent the uniaxial magnetic 
anisotropy of the molecule {D > and B > 0). 
The parameters, D and B, have been determined by 
high frequency EPR and neutron scattering experiments 
[7|, lia 13]. This spin-Hamiltonian, in a semiclassical 
view, describes a double potential well, in which oppo- 
site projections of the spin onto the z-axis are separated 
by an anisotropy energy barrier ^ DS^ + BS'^ 60 K) 
(see fig. 1). The third term is the Zeeman energy due to 
the longitudinal component of an external magnetic field, 
Hz = Hcos9, where 9 is the angle between the field and 
the easy axis of the molecule. A magnetic field applied 
along the z-axis shifts the double well potential and the 
energies of the projections of the magnetization. At cer- 
tain values of the 2;-axis field (resonance fields) the levels 
TO and to' with antiparallel projections onto the z-axis 
have nearly the same energy, Hk ~ kD/g^B = 0.44 T 
(^k = m + to') (see fig. 1). 

At these resonances, MQT is turned on by interac- 
tions that break the axial symmetry and mix the levels to 
and to'. These off-diagonal terms have different origins: 
(a) Ht is the Zeeman interaction associated with the 
transverse component of the external magnetic field; (b) 
transverse anisotropy terms are included in 7i_4; and (c) 
H' characterizes the transverse components of magnetic 
fields due to inter-molecular dipolar interactions and hy- 
perfinc nuclear fields. Off-diagonal terms in the Hamil- 
tonian lift the degeneracy of the spin-levels at the res- 
onances creating an energy difference between symmet- 
ric and antisymmetric superpositions of spin-projections 
that is known as the tunnel splitting, A (see fig. 1). 

Figure 2 shows magnetization curves measured at 0.6 
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FIG. 3: 3D representations of the anisotropy barrier. (A) 
Uniaxial anisotropy barrier in the absence of transverse 
anisotropy terms. The colored arrows represent the preferred 
orientation of the spin along the z-axis. (B) Anisotropy barri- 
ers corresponding to opposite signs of second order anisotropy, 
±E{S^ — Sy). The lines represent the hard (H) and medium 
(M) axes which are separated by 90 degrees. A change in the 
sign of E corresponds to a rotation of the transverse magnetic 
ax;es by 90 degrees. (C) Anisotropy barrier due to a fourth 
order anisotropy, (7(5* +5*1). In this case, there are two hard 
and two medium transverse axes separated by 45 degrees. A 
change of the sign of C produces a 45 degree rotation of the 
transverse magnetic axes. 

K for several sweep rates of the longitudinal field, a — 
dH/dt, for a single crystal of Mni2-ac. The abrupt in- 
creases of the magnetization toward the equilibrium mag- 
netization {M/Ms = 1) are due to MQT at the resonant 
fields. It is important to note for our discussion that 
MQT is observed at consecutive resonances. This has 
important implications for the understanding of MQT, 
because transverse anisotropy terms introduce selection 
rules and the only interaction that allows MQT at odd- 
resonances (k = 1, 3, 5...) is a transverse magnetic field. 

A. Transverse Interactions 

The tunneling characteristics depend on the form of 
the off-diagonal terms. In this subsection we will examine 
the consequences of the off-diagonal terms in the Hamil- 
tonian on MQT. We consider both transverse magnetic 
field and anisotropy terms. 

Magnetic Fields. The simplest expression for the off 
diagonal part of the Hamiltonian of Eq. (1), involving a 
transverse magnetic field, Ht, is 

Ht = —9^J■BHT{SxCos(f) + SySincf)) . (2) 

This represents the Zeeman energy for a field in the 
x-y plane at an angle (j) with respect to the x-axis. 
The tunnel splitting. A, is very sensitive to this field, 

Ht = \J H"^ + H"^ ■ The dependence of A on Ht for small 

transverse magnetic fields [Ht ^ Ht, = 2DS/giiB, the 
anisotropy field) is given by )32j. 

Ak{HT) ^ gk ( (3) 
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FIG. 4: Tunnel splitting of the ground state resonance, k — 0, 
as a function of the orientation (0) of a transverse field of 0.3 
T applied in the hard magnetic plane of a molecule for (A) 
only second, and (B) only fourth order anisotropy terms. We 
have used D = 548 mK, B = 1.17 mK, £ = 10 mK and C = 
0.022 mK in Eq. (1). The graphic clearly shows the 2-fold 
and 4-fold symmetries imposed by these different anisotropy 
terms. 

where gk = [(2S~k~i)]^ ^ y^ (2S-fc)!(2g )l = 25 - k. 

The power law dependence of A on Ht causes the tunnel 
splitting to vary by many orders of magnitude for trans- 
verse fields in the range of a few Tesla. This allows the 
study of MQT with a wide range of experimental tech- 
niques that go from quasi-static magnetization measure- 
ments (A//i ~ Hz) to high frequency EPR experiments 
{A/h 100 GHz) simply by varying the magnitude of 
the applied transverse field. 

Transverse Anisotropies. Now we examine the effect 
of transverse anisotropy terms in the Hamiltonian. Con- 
sidering only the lowest order terms, 7Y_4 has the form 

Ha^E{SI-SI) + C{SX+S^_) . (4) 

The first term is the second order anisotropy which is 
allowed for SMMs with rhombic symmetry (e.g. Feg [3^ 
and, as will be shown below, is also present in Mni2- 
ac due to disorder that lowers the S4 symmetry of the 
molecules 16. 18, 19 1 ) . The second term is fourth order in 
the spin operators and is the lowest order anisotropy al- 
lowed with tetragonal symmetry. This form of the trans- 
verse anisotropy has been observed in EPR studies of 
several SMMs of tetragonal symmetry such as, Mni2-ac 
[H, Mni2-BrAc and Ni4 [3. 

Let us consider in detail the consequences of each 
transverse anisotropy term. In fig. 3 we show the shape 
of the classical anisotropy barrier separating antiparallel 
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orientations of the spin z-projections (colored arrows). In 
the absence of transverse terms (fig. 3A) , the anisotropy 
barrier is determined by the uniaxial anisotropy of the 
molecules (first and second terms in Hamiltonian of Eq. 
(1)). These uniaxial anisotropy terms determine a hard 
anisotropy plane between opposite orientations of the 
spin states along the easy magnetic z-axis. Note that 
this barrier is isotropic in the x-y plane. In this case, 
the tunnel splitting does not depend on the orientation 
of the transverse field in the hard plane. 

This situation changes in the presence of a transverse 
anisotropy. Fig. 3B shows how the anisotropy barrier is 
modified by a second order anisotropy term of the form 
E{Sl. — Sy). This term introduces one hard and one 
medium axis in the hard plane (the x-y plane) that are 
separated by 90 degrees. That is, for positive E the x is 
hard and y axis is medium. A change of the sign of E 
leads to a 90 degree rotation of these axes. In this case, 
the tunnel splitting. A, depends on the azimuthal an- 
gle, 0, of the apphed transverse field, Ht- A transverse 
magnetic field applied along the medium axis produces a 
larger tunnel splitting than the same field applied along 
the hard axis. This leads to an oscillatory dependence of 
A on with 2 maxima and minima separated by 90 de- 
grees (see fig. 4A). A second order anisotropy also intro- 
duces MQT selection rules. In the absence of a transverse 
field, this term only allows MQT for resonances that are 
even (i.e., k — 2i with i an integer). 

Fig. 3C shows the modification of the anisotropy bar- 
rier in the presence of a fourth order anisotropy term 
of the form C(5^]l + 5*1). This anisotropy produces two 
medium and two hard axes in the hard plane. The sep- 
aration between medium and hard axes is 45 degrees. 
Thus the tunnel splitting will have a 4-fold pattern of os- 
cillations as a function of the angle of a transverse field, 
with maxima and minima spaced by a 45 degrees (see 
fig. 4B). This term introduces a selection rule that al- 
lows tunneling transitions from the ground state for res- 
onances that are a multiple of 4 (fc = Ai). In fig. 4 we 
show the expected behavior of the tunnel splitting of the 
ground state for resonance fc = versus the angle of the 
applied transverse field, Ht — 0.3 T. The tunnel splitting 
is calculated using the Hamiltonian of Eq. (1) with D = 
548 mK and B — 1.17 mK and considering only second 
order anisotropy {E = 10 mK, fig. 4A) and only fourth 
order anisotropy (C = 0.022 mK, fig. 4B). 



B. Disorder 

Recent experimental results have shown that MQT 
in Mni2-ac is modulated by off-diagonal terms that are 
generated by disorder. Disorder allows for anisotropy 
terms that are lower order in the spin-operators than 
those imposed by the average molecule site symme- 
try in SMM crystals. Disorder can also lead to tilts 
of the eas y a xes from the crystallographic easy axis 
H m 111 "El E in 111. The phenomena that have 



been explained by disorder can be summarized as follows: 
a) observation of MQT relaxation at fc-resonances that 
are not allowed by the quantum selection rules imposed 
by the symmetry of the molecule and b) non-exponential 
magnetic relaxation which suggests the existence of a dis- 
tribution of timnel splittings. 

Two distinct models of disorder have been proposed. 
Chudnovsky and Garanin jl^ il5| proposed that random 
line dislocations in a crystal lead, via magnetoelastic in- 
teractions, to a lower molecule symmetry and a broad 
distribution of tunneling rates. Subsequent magnetic re- 
laxation experiments indeed showed the existence of a 
broad distribution of tunneling rates and were analyzed 
in terms of this model 9, 10]. In contrast, Cornia et 
al. suggested, based on detailed x-ray analysis, that 
variations in the position of the two hydrogen-bonded 
acetic acid molecules surrounding the Mni2-ac clusters 
lead to a discrete set of molecules with lower symmetry 
than tetragonal. More recent magnetic relaxation experi- 
ments [H and high frequency EPR experiments |3 have 
confirmed the latter model, showing a 2-fold symmetry 
of the tunnel splitting as a function of the direction of an 
external transverse field. Moreover, the observation of a 
distribution of transverse fields in Mni2-BrAc |2| and 
tilts in Mni2-ac ,27J suggests that tilts of the easy axes 
of the molecules are responsible for the MQT relaxation 
at odd-A; resonances; in the case of Mni2-ac, these tilts 
are caused by the solvent disorder. We summarize these 
models and their consequences below. 

Line Dislocations. Chudnovsky and Garanin '14', TB'l 
considered a random distribution of line dislocations with 
collinear axes to calculate a representative distribution 
of second order transverse anisotropics. They found: 
a) the corresponding distribution of tunnel splittings is 
broad on a logarithmic scale; b) the mode of the dis- 
tribution is at -E = 0; and c), for small concentrations 
of dislocations, most molecules are far from the dislo- 
cation cores and the magnetic axes associated with the 
transverse anisotropy are oriented randomly. The first 
point can explain the non-exponential relaxation found 
in Landau-Zener relaxation experiments carried out at 
several MQT resonances, and at different longitudinal 
magnetic field sweep rates . Furthermore, experiments 
in which a Mni2-ac single crystal was treated with rapid 
thermal changes could be explained in terms of the dis- 
location model Point (b) leads to a distribution of 
tunnel splittings, with a very long tail towards small val- 
ues of the tunnel splitting. In a more recent experimental 
study, a Landau-Zener method that involved crossing the 
same MQT resonance several times permitted the study 
of nearly the whole distribution of tunnel splittings in a 
Mni2-ac single crystal at a single resonance J^. The 
results obtained were compared to the distribution func- 
tion proposed by Chudnovsky and Garanin. A distribu- 
tion of second order anisotropics with mode at _E = was 
not able to explain the relaxation data. More recently, 
both magnetic relaxation and high frequency EPR ex- 
periments 13 confirmed the discrete nature of the 
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distribution of transverse anisotropies and, moreover, the 
fact that the magnetic axes of the transverse anisotropy 
is also discretely distributed along particular directions 
in the hard anisotropy plane of the molecules, contrary 
to the random distribution expected from the dislocation 
model. 



Solvent Disorder. From x-ray diffraction studies Cor- 
nia et al. showed that the fourfold symmetry of 

Mni2-ac molecules is lowered by disorder of the acetic 
acid molecules (2 molecules with 4 possible sites). This 
gives rise to a set of six different molecules (4 of them 
with E 0), with 7/8 or 88% of molecules thus having 
a non-zero second order transverse anisotropy. The fact 
that the second order anisotropy is generated by disorder 
in the solvent molecules implies that both the abundance 
and magnitude of the E for each isomer could depend 
on the synthesis process as well as on the solvent losses 
that a particular crystal experienced prior to measure- 
ment. From x-ray data Cornia et al. calculated the E 
parameters for the different isomers These should 

be taken as estimates, since the data were obtained from 
one crystal and the E-parameters were computed using 
an empirical model. Nonetheless, this model represents 
an important step in the understanding of MQT in Mni2- 
ac. 



The relevant points of the solvent disorder model can 
be summarized as follows: a) there is a discrete set of 
E values in a sample. This means that the distribu- 
tion of tunnel splittings is a discrete function with peaks 
at several E values, b) There are equal popula- 
tions of isomers having opposite signs of the second or- 
der anisotropy ±E. This means that the medium and 
hard axes associated with each E value will be discretely 
distributed over the hard plane of the crystal with a sep- 
aration of 90 degrees. And c) the hard axes of the second 
order anisotropy are rotated away the hard/medium axes 
corresponding to the fourth order anisotropy; this rota- 
tion angle depends on the precise details of the solvent 
structure, and is ~30 degrees for Mni2-ac.. The first 
point has been confirmed by high frequency EPR experi- 
ments by the observation of discrete peaks in the absorp- 
tion spectra that correspond to different discrete E val- 
ues and have the angular dependence expected from this 
model [isj . The second point has been confirmed through 
magnetic relaxation measurements in which a portion of 
molecules having different signs of E have been studied 
[TtI . The third point has been observed in EPR experi- 
ments and has implications for the MQT discussed below. 
We note that both models suggest the presence of tilts 
of the easy magnetic axes of the molecules. These latter 
results, which have recently been confirmed via density 
functional calculations , will be presented in detail in 
the experimental sections of this paper. 



C. The Effect of Disorder on MQT 

In this section we will present model calculations of 
the dependence of the tunnel splitting on the angle and 
the magnitude of an external transverse field taking into 
account the presence of both second and fourth order 
anisotropies. These terms must be considered on an 
equal footing because they produce comparable tunnel 
splittings. This is seen as follows. In perturbation 
theory, the splitting between levels m and m' associ- 
ated with the second order anisotropy is approximately 
Ae D{E/2DY"''^"''^^^, while that associated with 
the fourth order anisotropy is Ac ~ 0(0/20)''-^'-"''^/'^. 
Thus, provided E^ ^ DC these are comparable in mag- 
nitude, independent of m and m' . This is the case for 
Mni2-ac, asE ^ IQ-^ K, C ~ IQ-^ K and D ~ 1 K. For 
this reason there is an intricate interplay between these 
anisotropy terms that also occurs over a broad range of 
transverse magnetic fields, and can thus be probed both 
in EPR and MQT studies. 

In Mni2-ac these 2nd and 4th order transverse 
anisotropy terms have a different origin. The fourth order 
transverse anisotropy term is imposed by the symmetry 
of the molecule, while the second order term is generated 
by a distribution of solvent micro-environments. How- 
ever, there is no reason to assume that these anisotropies 
are commensurate since they come from different sources; 
one is associated with the ideal S4 symmetry, the other 
is disorder induced and lowers this symmetry. Therefore, 
we will consider a misalignment angle, /3, between the 
magnetic axes associated with each transverse anisotropy 
term. 

We write the transverse anisotropy part, 7i^, of the 
Hamiltonian (Eq. (4)) as follows: 

E{ei{Sl-Sl) + e2iS,-Sy + SyS,)) + C{Sl+St) , (5) 
where, 

ei = {cos^ (3 — sin^jS) 62 = 2cosf]sin(3 . (6) 

This introduces an angle (3 between hard axes of the sec- 
ond and fourth order anisotropy terms (see fig. 5), the 
'i?' and 'C terms respectively. 

In the next subsections we will show the effect of hav- 
ing a misalignment (/3) on the symmetry of MQT and 
on the transverse field dependence of the tunnel splitting 
for different orientations of the applied field. The cal- 
culations have been done for parameters and conditions 
close to those shown in the experimental sections of this 
article. In the following we consider in some detail the 
consequences of the solvent disorder model as this can ad- 
equately explain many of our experimental observations. 

As mentioned, one of the important consequences of 
the solvent disorder model is the fact that there must be 
equal populations of molecules having opposite signs of 
E. This is because, on average, Mni2-ac crystals have S4 
or tetragonal site symmetry. Note that fig. 5 represents 
the picture for the molecules having E > 0. Thus there 
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FIG. 5: 3-D representation of the anisotropy barriers corre- 
sponding to fourth (upper graph) and second order (lower 
graph) anisotropies with an angle /3 between the correspond- 
ing hard axes. 
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FIG. 6: Behavior of the ground state splitting, Afc=o, versus 
the orientation of a transverse field, Ht = 9 T, for /3 = and 
different values of i? > 0. 

will be molecules in the sample with the magnetic axes of 
E rotated by 90 degrees with respect to those represented 
in fig. 5. 

Commensurate Transverse Interactions. We consider 
the case of i? > and /? = in the presence of large 
transverse field, as this field corresponds to that used in 
high frequency EPR experiments (section IV). Fig. 6 
shows the behavior of the ground state tunnel splitting, 
^10,-10 (resonance k = Q), versus the angle, 4>, of an 
external transverse field of 9 Tesla. The calculations have 
been done for /3 = 0, £> = 548 mK, B = 1.17 mK, 
C = 0.022 mK and different values of E. 

One interesting result in this figure is the observation of 
fourfold pattern of maxima imposed by the fourth order 
anisotropy. There are four maxima in A that occur very 
close to the medium axes of C for small values of E. Only 
when E becomes very large do the maxima move towards 



£=3mK ---- E=l-3mK 
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FIG. 7: (a) Ground state tunnel splitting, Afc=0) versus the 
orientation of a transverse field, Ht = 9 T, for /3 = and 
E = -1-3 mK (solid line) and £; = -3 mK (dashed line), (b) 
Average value of the tunnel splitting assuming both iJ-signs. 
(c) Difference between the tunnel splitting values of _E > 
and E <Q. 

the medium axis of E to give, at high enough i?-values 
(i? >> 20 mK), a twofold pattern of maxima. However, 
for the parameters estimated with the solvent disorder 
model, ii^ < 10 — 15 mK, the tunnel splitting is expected 
to show fourfold maxima in the direction of the medium 
axes of C. 

Fig. 7a shows the behavior of the ground tunnel split- 
ting, Afc=0) for the same conditions as fig. 6, that is for 
/? = 0, but taking into account the equal population of 
molecules having opposite signs of E. In this case we 
have used E = Eav = ±3 mK. The modulation of the 
fourfold symmetry by E is opposite for different signs of 
E. However, the maxima are approximately at the same 
position, independent of the sign of E. In fig. 7b we 
show the (t> dependence of the average value of the tun- 
nel splitting (A£;>o -I- Ae<o)/2. This average exhibits a 
symmetric four-fold rotation pattern, the same as would 
be generated via only a fourth-order anisotropy. There- 
fore, a measurement of the average tunnel splitting of a 
crystal cannot distinguish between second and fourth or- 
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FIG. 8: (a) Ground state tunnel splitting, Afc=0) versus the 
orientation of a transverse field, Ht = 9 T, for E = +3 mK 
and different values of /3. (b) Average value of the tunnel 
splitting assuming both i?-signs for each value of /3; all of 
the curves lie on top of each other, (c) Difference between 
the tunnel splitting values for E > and E < for different 
/3- values. 

der transverse anisotropies if there are equal populations 
of molecules with opposite signs of E. 

The absolute difference between the tunnel splittings, 
corresponding to opposite signs of E, is shown in fig. 7c. 
In this case, (3 = 0, the difference is maximal along the 
hard axes of C and vanishes along the medium axes of C. 
Importantly, this provides a means of inferring the pres- 
ence of a second-order anisotropy through an appropriate 
experiment (this will be shown below and in section IV) . 

Incommensurate Transverse Interactions. We now 
consider the effect of a misalignment, /3 7^ 0. Fig. 8a 
shows the behavior of the ground tunnel splitting, A^^o > 
versus the angle of the applied transverse field, Ht = 9 T 
for E = 3 mK and different values of /3, from f3 = —45° to 
P = 45° in increments of 15 degrees. Note that the black 
line, /3 = 0, is the result presented in fig. 7a. A misalign- 
ment P ^ generates an asymmetry between the maxima 
of the tunnel splitting. For example, for P = —45°, the 
maximum at = 45° is bigger than the maximum at 
(j) = 135°. Note that the hard axis of E (HE) for this 



value oi P is (j) = —45° -|- nl80° (with n an integer) while 
the medium axis (ME) is along </) = 45° -I- nl80°. In gen- 
eral, even though the ME axis is at = /3 + (2ri + 1)90°, 
the maxima of the tunnel splitting are in the direction of 
the medium axes of C . The second order anisotropy in- 
troduces an asymmetric modulation of the maxima and 
minima of the tunnel splittings that depends on the value 
and the sign of E and on the misalignment angle p. 

Fig. 8b shows that the average value of the tunnel 
splitting is independent of the angle of misalignment P 
and has four maxima. All the results collapse in the 
same curve. So again a direct measurement of the av- 
erage value of A would not give information about the 
misalignment between these anisotropies. 

Fig. 8c shows the difference between A£;>o and /S.e<o 
for the same parameters used in the previous calculation 
and different values of P from —45° to 45°. This differ- 
ence has four maxima for all the angles p. However, the 
positions of the maxima are different for different P val- 
ues. In fact, the position of the maxima depends directly 
on the value of P as (j)max = P + n90°. Consequently, 
a measurement of this difference would not only provide 
the value of E but also give the angle of misalignment 
between second and fourth order transverse anisotropies. 

We thus note that within the solvent disorder model 
the tunnel splitting of resonance fc = (the ground 
state degeneracy) at high transverse field values will 
show a fourfold rotation pattern with respect the angle 
of orientation of the applied field. This symmetry is 
modulated by the second order anisotropy. However, 
equal populations of molecules with opposite signs of E 
will give an average tunnel splitting with four maxima 
that is indistinguishable from that corresponding fourth 
order transverse anisotropy. The only way to determine 
the presence of second order anisotropy is by examining 
the difference of the tunnel splittings corresponding to 
molecules with opposite signs of E. The latter would 
give the value of E and the misalignment angle between 
the second and fourth order anisotropies. We will show 
how such measurements are possible via two different 
techniques in the experimental part of this paper. 

Incommensurate Transverse Interactions and MQT. 
Let us now consider the situation relevant to understand- 
ing the physics of MQT in magnetic studies in which 
much smaller tunnel splittings are probed 10~^ K). 
We will analyze the behavior of the splitting for reso- 
nance fc = 6 (m = — 10,m' = 4) versus the angle, (/), 
of a small external transverse field. The calculation has 
been done with the same D, B and C parameters of the 
Hamiltonian as those of the previous calculations. We 
have used jfi'l = Eav = 3 mK and a transverse magnetic 
field, Ht = 0.35 T, which corresponds to the situation 
studied in the dc Landau-Zener relaxation experiments 
that will be presented in section III. 

The results are shown in fig. 9. This situation is sub- 
stantially different from that of the previous case of res- 
onance fc = 0. The tunnel splitting has a pattern that 
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FIG. 9: Ground state tunnel splitting, Afc^g, versus the ori- 
entation of a transverse field, Ht = 0.35 T, for E = +3 mK 
and different values of /3. 




^ 45 90 135 180 

(z> (degrees) 

FIG. 10: Average value of the ground state tunnel splitting, 
Afc^g, for opposite signs of E {\E — 3| mK), versus the ori- 
entation of a transverse field, Ht = 0.35 T, and for different 
values of p. 



goes from two-fold maxima for /? = ±45° to four-fold 
maxima for /3 = 0. Interestingly, the maxima are in 
the directions of the medium axes of C independent of 
the direction of the medium axis of E. For example, for 
(3 = —30° the two maxima of the tunnel splitting are 
at (t>max = 45°, 135° while the directions of the medium 
axis of are 0M£; = /3 + 90° = 60°, 240°. So, for some 
/3-values, the tunnel splitting exhibits a twofold pattern 
of maxima with position determined by the fourth order 
anisotropy. We have done these calculations for bigger 
values of E (not shown in this paper) that indicate that 
the range of /3-values around /3 = that exhibits fourfold 
symmetry is narrower the bigger the value of E. However, 
for i?-values smaller than 30 mK the maxima positions 
are still determined by the fourth order anisotropy. This 
constitutes an unexpected result in a system with second 
and fourth order anisotropics that is first pointed out in 
this work and has important consequences for the inter- 
pretation of the experimental measurements that have 
been previously published by some of the authors of this 
work |TM IT9I I. as will be explained in the experimental 
sections III and IV. 

We have calculated the average value of the tunnel 
splitting for different angles (3 assuming an equal popu- 
lation of molecules with different signs of E. The results 
are shown in fig. 10. For all values of (3 there is clearly a 
fourfold rotation pattern of maxima of the average tun- 
nel splitting with symmetric maxima and minima along 
directions determined by the fourth order anisotropy. 
Again, the only method to determine the E value and 
the relative orientation between C and i?, /3, is to study 
a subset of molecules having only one sign of E. In the 
latter case, depending on the value of (3, there is a possi- 
bility that the molecules selected will also show fourfold 
symmetry (e.g. if /3 ~ 0). It would then not be possi- 
ble to conclude that such molecules have a second order 
anisotropy. As we will show in sections III and IV this is 
not the case experimentally. We will show that a selection 



of a subset of molecules with one sign of E shows twofold 
symmetry which indicates that the angle (3 of misalign- 
ment is close to (3 ^ ±45° (specifically, (3 — —30°). 

Disorder and Berry Phase Effects. Quantum phase 
interference is one of the most important phenomena ob- 
served in SMMs. Interference effects (Berry phase) in 
MQT where first discovered by Loss [33 and calculated 
for a nanomagnet with a biaxial anisotropy by Garg |37j . 
This phenomenon is due to interference of the quantum 
tunneling trajectories of the magnetization and has been 
clearly observed in two SMMs to date 33, 38] . The first 
observation of the Berry phase was by Wernsdorfer and 
Sessoli (sil in the Fes SMM. Fcs has both second and 
fourth order transverse anisotropy. The observation of 
quantum oscillations was done by applying a transverse 
field along the direction of the hard axis of the second or- 
der anisotropy (HE), which in this case also corresponds 
to the direction of one of the hard axes of the fourth 
anisotropy term, (3 — Q. Fig. 11 shows the dependence 
of the tunnel splitting of resonance A: = on the external 
field applied along cf) = Q (hard E anisotropy axis) and 
<j) — 90° (medium E anisotropy axis). We have used the 
parameters given in reference |33l |. 

The field spacing between quantum tunneling oscilla- 
tions of fig. 11 can be described in terms of the anisotropy 
parameters D and E of the Hamiltonian using a semiclas- 
sical approach 37], 



AH^—^2E{E + D) , (7) 

This is ^ 0.23 T for Fes which is smaller than the spacing 
resulting from the exact diagonalization of the Hamilto- 
nian including fourth order transverse anisotropy, AH — 
0.41 T, and was observed in reference [s^. The reason 
is that eq. (7) only considers the presence of a second 
order anisotropy (23 ■ Note that the authors of ref. 
assumed collinear anisotropics to fit the data, thus the 
angle between the HE and HCl axes is /? = (see fig. 




FIG. 11: Feg ground state tunnel splitting, A^^o, versus 
the magnitude of a transverse field applied along the hard 
(HE) and the medium (ME) anisotropy axes. The parame- 
ters used in this calculation were taken from reference |3^. 
The drawings above represent the x — y plane projections of 
the anisotropy barrier for Feg in the presence of a transverse 
field applied at different angles 0. The white lines show two 
hypothetical quantum tunneling trajectories. When Ht is 
applied along the hard anisotropy axis, = 0, the barrier 
remains symmetric with respect to the field. In this case, the 
trajectories interfere. For transverse fields not aligned with 
the hard axis, an asymmetric distortion of the barrier leads to 
non-equivalent MQT trajectories, destroying the interference. 

5). We also have used coUinear second and fourth or- 
der anisotropy terms in the calculated data shown in fig. 
11. So we can see that the effect of having a fourth 
order anisotropy that is coUinear with a second order 
anisotropy, /3 = 0, only modifies the pattern of oscilla- 
tions but not their structure and shape. Quantum tun- 
neling oscillations were also observed in [Mni2]~^ which 
only has a second order transverse anisotropy [s^ . In 
this case, the spacing between oscillations was given by 
eq. 7. 

Park and Garg [40l | calculated the quantum tunneling 
oscillations in a system with only fourth order transverse 
anisotropy using the Hamiltonian of Mni2-ac (Eq. (1)). 
Fig. 12 shows calculations of the tunnel splitting for 
different resonances by using eq. (1) with D — 556 mK, 
B = 1.1 mK and C = 0.03 mK. 

In order to take into account the effect of a mis- 
alignment, /3 7^ 0, between second and fourth order 
anisotropies we have calculated the dependence of the 
ground state tunnel splitting on the magnitude of an ex- 
ternal transverse field applied along different character- 
istic directions in the hard plane of a molecule. For this, 
we have used (3 = -30°, D = 548 mK, B = 1.17 mK and 
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FIG. 12: Transverse field dependence of the tunnel splittings 
for Mni2-ac , for several resonances, assuming only fourth 
order anisotropy. The transverse field is applied along one of 
the hard axes of C (HC). The drawings above the figure show 
the distortion of the anisotropy barrier due to a transverse 
field applied at different angles with respect to the hard axes 
of the fourth order anisotropy. White lines represent different 
tunneling trajectories. 

C = 0.022 mK and different values of i? > 0. For clarity, 
this situation corresponds to having the following direc- 
tions for the characteristic transverse anisotropy axes: 
(t>HC = 0,90°, 180°, 270°, (t>Mc = 45°, 135°, 215°, 305°, 
(l)HE = /3 + nl80° = -30°, 150° and = 60°, 240°. 

Examining the result for (3 — —30° in fig. 9 (red 
curve) one can see that the tunnel splitting has twofold 
symmetry with maxima at </> = 45°, 225° and minima 
at = 135°, 315°. These directions correspond to the 
medium axes of C (MCI and MC2). We have thus calcu- 
lated the dependence of the ground state tunnel splitting 
for resonances fc = 5, 6 and 7 for the field applied along 
MCI and MC2. These are shown in fig. 13a, 13b and 13c 
respectively, for a transverse field applied along <j) — 45° 
(tunnel splitting maximum, MCI axis) and (f> = 135° 
(tunnel splitting minimum, MC2 axis) using E — 10 mK 
(thin lines) and £' = 15 mK (thick lines). Fig. 13 shows 
how the minimum in the tunnel splitting moves to higher 
fields and becomes deeper as E increases. 

The results show a very different structure and shape 
of the Berry phase oscillations as compared to the case 
when both anisotropies are commensurate. Note that the 
orientations we used for the applied transverse field corre- 
spond to the medium axes of the fourth order anisotropy 
where one does not expect to have Berry phase oscilla- 
tions. Moreover, these two orientations, </> = 45° and 
(f) = 135°, do not coincide with the hard axis of the 
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FIG. 13: Transverse field dependence of tlie ground state tun- 
nel splitting of resonances k — 5, k — 6 and k = 7 for Mni2- 
ac, with E = 10 mK and E = 15 mK. The transverse field 
is applied along the = 45° and — 135° directions, which 
correspond to the positions of the maxima and minima of A 
for a misalignment angle /3 = —30°. 

second order anisotropy, which for f3 = —30° is along 
(j) = -30° + nl80° . 

In order to have a more complete picture of the effect 
of second and fourth anisotropics with non-colhnear axes 
on the Berry phase phenomena we have calculated the 
dependence of the ground state tunnel splitting of reso- 
nance fc = 7 on the magnitude of a transverse field ap- 
plied at different angles, from = 45° (which correspond 
to one of the medium axes of C) to (/) — 225° (which 
correspond to the opposite orientation of the field along 
the same hard C-axis). The results are shown in a color 
contour plot in fig. 14. The first thing to point out is 
that the tunnel splitting still has zeros in this situation. 
However, the most significant fact is that these zeros do 
not appear at an angle characteristic of the transverse 
anisotropics. Moreover, the structure, shape and posi- 
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FIG. 14: Color contour plot of the transverse field dependence 
of the ground state tunnel splitting of resonance k = 7 for 
different angles (p for a misalignment, (3 — —30°, between the 
hard anisotropy axes of E and C. The vertical lines represent 
the orientations of hard and medium axes of the second and 
fourth order anisotropy terms. This misalignment generates 
a new and interesting pattern of Berry phase zeros that does 
not correspond to any of the anisotropics separately. 




FIG. 15: Projection of the anisotropy barrier onto the x — y 
plane for second (right-upper graphic) and fourth (left-upper 
graphic) order transverse anisotropy terms. The addition 
of both anisotropics leads to an asymmetric barrier (lower 
graphic) that makes it difficult to see the orientations of a 
transverse field that would generate equivalent quantum tun- 
neling trajectories and, therefore. Berry phase interference 
phenomenon. 

tion of the zeros are completely independent of the Berry 
phase corresponding to each anisotropy term separately. 
Note that the calculations shown in fig. 13 correspond to 
transverse fields applied along (f> = 45° and 135° which 
correspond to to the directions of the maximum and min- 
imum values of the tunnel splitting of fig. 9 respectively. 
In the case oi (f) — 135° the tunnel splitting is close to 
one of the zeros (</) ~ 125° and Ht ^ 0.4 T in fig. 14) 
but far from all the others, explaining the observation of 
only one incomplete oscillation in fig. 13a. This new and 
unusual structure of the Berry phase zeros can be better 
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understood by looking at the graphic representation of 
the anisotropy barrier of fig. 15. The addition of sec- 
ond (right- upper illustration in fig. 15) and fourth order 
(left-upper) transverse anisotropics in a SMM leads to an 
asymmetric barrier (center-lower) where symmetry does 
not permit a direct identification of the field that gener- 
ates equivalent quantum tunneling trajectories that can 
interfere. 

From the results, we can conclude that the combina- 
tion of incommensurate transverse anisotropy terms of 
different order in the spin-operators can lead to an inter- 
esting situation in which the resulting magnetic response 
does not depend in any simply way on the form of either 
anisotropy term separately. The parameters used in the 
above simulation were chosen because they are within the 
range of values that can explain the experimental results 
presented in sections III and IV. 



III. MAGNETIC RELAXATION 
EXPERIMENTS 

We have carried out magnetic relaxation measure- 
ments in a single crystal of deuterated Mni2-ac in the 
pure quantum regime (T = 0.6 K) in which relaxation is 
by MQT without thermal activation over the anisotropy 
barrier |^. Deuterated crystals were studied because 
the purity of the chemicals used in the synthesis leads to 
very high quality crystals . 

We have used a high sensitivity micro-Hall effect mag- 
netometer to measure the magnetic response of a 
Mni2-ac single crystal of ~ 100 micrometer size and nee- 
dle shape. We measure the longitudinal component of 
the magnetization of the sample (z-component) by plac- 
ing the crystal with one of its faces parallel to the sensor 
plane and one end just over the cross point of the micro- 
Hall sensor. The magnetometer was placed inside a low 
temperature He'^ system. A superconducting vector-field 
magnet was used to apply high magnetic fields at arbi- 
trary directions with respect to the crystal axis. 



A. Landau-Zener Method 

The Landau-Zener (LZ) method has been used to 
study quantum tunnel splittings in SMMs 's^ and has 
become a powerful tool to check for distributions of 
dipolar and nuclear interactions [s^, molecular micro- 
environments_P, ITol llTl Il8|| or internal transverse mag- 
netic fields 1221 in these materials. The method consists 
in crossing a MQT resonance by sweeping the longitudi- 
nal magnetic field at a constant rate, a — dH/dt, and 
measuring the fractional change of the magnetization in 
the process. The anti-crossing of the spin levels m and m' 
of resonance k = m + m' is shown in the inset of fig. 16. 
For an ideal system of non-interacting and monodisperse 
SMMs, and for low enough temperatures (where thermal 
relaxation is negligible), the normalized change of magne- 




H(T) 

FIG. 16: Landau-Zener multi-crossing experiment in a Mni2- 
ac single crystal measured at 0.6 K by sweeping the longitu- 
dinal magnetic field at a constant rate, a — 6.6 x 10"'' T/s, 
multiple times across the k = 7 resonance. The inset shows 
a representation of the energy levels, m and m', at the anti- 
crossing point. 

tization, {Mbefore - Mafter)/iMbefore - Meq), is related 
to the probability for a molecule to reverse its magne- 
tization by quantum tunneling. The bigger the MQT 
probability the larger the magnetization change will be. 
This MQT probability is related to the tunnel splitting, 
A, by the LZ formula |4||, 

PLZ = l-exp[-- , 8 

where vq = gj-isi'^S — k) and v^a is the energy sweep 
rate. Rlz = 1 ^ Plz is the probability for a molecule 
to remain in the mctastable state | m') after crossing the 
resonance. 

It is important to note that this relation is only valid 
if the internal energy sweep rate a single molecule expe- 
riences is proportional to the external sweep rate of the 
magnetic field. This is not satisfied if there are inter- 
nal dipolar or nuclear fields. In fact, it has been shown 
in Fes |H| and Mni2-ac [l3l that changing dipolar fields 
lead to deviations from the LZ formula. In order to avoid 
these effects the external sweep rate must be fast enough 
to have small magnetization change (i.e. small changes 
of dipolar fields) in the crossing process. The critical 
lower value of the magnetic field sweep rate, ac, that is 
needed to avoid this situation has been determined to be 
'p^s foj- Mni2-ac JTJ- Due to this, all the experi- 
ments presented in this section have been conducted with 
a > ac- 

B. Multi-crossing Landau-Zener Measurements 

When there is a distribution of quantum splittings in 
the sample each molecule has a different MQT proba- 
bility and the relaxation of the magnetization should re- 
flect this fact. In this case, the MQT probability depends 
on the distribution of tunnel splittings of the molecules 
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FIG. 17: Landau-Zener probability to remain in the 
metastable well in multi-crossing relaxations of resonances 
= 6, 7 and 8 using different sweep rates (3.33 x 10"^ to 
1.33 X 10~^). All the measurements were performed begin- 
ning with a saturation magnetization at T — 0.6 K. 

that arc in the metastable well before crossing a reso- 
nance. After a crossing of a resonance, those molecules 
with the largest tunnel splitting values and, consequently, 
the highest MQT probability will represent the maximum 
contribution to the fractional change of the magnetiza- 
tion. Correspondingly, those molecules with smaller tun- 
nel splitting values will remain in the metastable well. 
Due to this, the LZ relaxation method can be used to 
determine the distribution of tunnel splittings in a sam- 
ple, and to select different parts of the distribution for 
independent study, as we will show in subsection III.C. 

In order to extract the complete distribution of tunnel 
splittings in Mni2-ac we have used a modification of the 
LZ method that consists of crossing a resonance multiple 
times, both for increasing and decreasing fields. As we 
have stated previously, molecules in the metastable well 
having the largest tunnel splitting values are most likely 
to relax in any given crossing of a resonance. Once a 
molecule has relaxed, it will no longer exhibit dynamics in 
subsequent crossings of the same resonance. Thus, only 
those molecules with the largest probability of tunneling, 
and which did not already tunnel, can contribute to the 
relaxation during subsequent crossings of the resonance. 
The repetition of this procedure many times enables a 
determination of the distribution of tunnel splittings in 
the sample over several orders of magnitude. We show 
an example of this multi-crossing LZ procedure for reso- 
nance fc = 7 in fig. 16. 

In a multi-crossing LZ measurement the probability to 
remain in the metastable well after crossing a resonance 
n-times is given by, 



Rlzv, = exp 



7rA2 1 



21^0 aeff 



(9) 



where cteff — a/n. If this expression describes the 
physics, then relaxation curves recorded at different 
sweep rates should scale when plotted as a function of 
the effective sweep rate. This can be clearly observed in 
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FIG. 18: Log-normal distribution functions extracted from 
the fit to the LZ relaxation curves of fig. 17. 

fig. 17, where we show LZ multi-crossing relaxation mea- 
surements of resonances k = 6, 7 and 8, carried out at 
different sweep rates a > ac (3.33 x 10"^ to 1.33 x 10~^). 
Small differences in the results were observed in three 
different crystals that were synthesized in the same way. 
These results clearly show that the MQT relaxation rate 
is not exponential and indicate the presence of a distribu- 
tion of tunnel splittings within the sample. These results 
confirm previous experimental observations of non expo- 
nential relaxation in Mni2-ac [ESS^S- Moreover, the 
large fraction of the relaxation that we are able to ob- 
serve with this method gives direct information on the 
width of the distribution of tunnel splittings. 

We have assumed a log-normal distribution of tunnel 
splittings to explain our observations. We take the form. 



fix) 



Aexp 



(10) 



where x = logA and Xc = logAc] Ac and a represent 
the center and the width of the distribution, respectively. 
The fit of the relaxation curves that are shown in fig. 17 
(solid lines) have been obtained by using. 



i?(a, n) 



Ruia,A,n)f(logA)dlogA 



(11) 



with Xc and a as free parameters. The resulting fits are 
in excellent accord with the experiments. The resultant 
distribution functions for resonances k — 6, 7 and 8 are 
shown in fig. 18. The center of the distribution increases 
with the resonance number, k, while the width remains 
almost constant for different resonances, being somewhat 
narrower for resonance k — 8. 

In section II. b we have discussed the two models that 
have been proposed to explain the presence of a broad 
distribution of tunnel splittings in terms of disorder. 
Before analyzing our results, note that a second order 
transverse anisotropy allows transitions for resonances 
2i, while a fourth order anisotropy, which is imposed by 
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ln{E/2D) 

FIG. 19: Distribution of the second order anisotropy parame- 
ter inferred from the log-normal distribution functions of res- 
onances k — 6 and 8 of fig. 18 using the expression of eq. 
(11). 

the symmetry of the molecules, only allows MQT relax- 
ation for resonance numbers which are a multiple of 4, 4z. 
Thus, a comparison between relaxation curves recorded 
at resonances fc = 6 and 8 should give us information 
about the origin of the tunnel splittings in this material. 
The relation between the tunnel splitting and the sec- 
ond order anisotropy is given by the next formula which 
follows from perturbation theory [33 |: 

ln(Afe/5,)/e, = hr(^^^ , (12) 

where gk and depend on k, S and D [s^ and were 
given in section II. a. Through this expression we can 
infer the distribution of second order anisotropy param- 
eters, f{ln{E/2D)), by taking the log-normal distribu- 
tions used to fit the relaxation curves of resonances k = 
6 and k = 8. The results are shown in fig. 19. The 
fact that both distributions do not scale when plotted 
as a function of ln{E/2D) indicates that second order 
anisotropy can not be the only origin of tunnel splittings 
within the sample. 

The distribution function of the second order 
anisotropy parameter predicted by the line dislocations 
model is given by [l5l |. 

1 / e^"' \ 
fhix) = —exp X = — , (13) 

where a; = ln£^ with E = E/2D. E^ is the width of the 
distribution of the anisotropy parameter E. E^ depends 
on the geometry of the crystal and on the concentration 
of dislocations per unit cell, c. Note that the mean value 
and width of the distribution given by this expression 
are not independent variables. By using this distribution 
with eq. (11) we have fit the relaxation curve recorded at 
resonance k — Q. The only free parameter in the fit is E^. 
The result is shown in fig. 20 (thin line) where the log- 
normal distribution extracted from our previous fit of the 
same relaxation curve has been included for comparison 
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FIG. 20: Tunnel splitting distribution of resonance k = 6 
expected from the line-dislocation model (thin line) compared 
with the log-normal distribution that fits the experimental 
data (thick line). 

(thick line) . We have chosen the mean value to be at the 
same position as that of the log-normal distribution. The 
value of Ec used to fit the data corresponds to a concen- 
tration of dislocations per unit cell of c 10~^. Clearly, 
the width of this distribution is many orders of magni- 
tude bigger than the log-normal distribution used to fit 
our data. This is due to the fact that the distribution of 
E predicted from the line dislocations model has a most 
probable value at i? = 0, which explains the long tail ob- 
served for low tunnel splitting values. Consequently, line 
dislocations would produce a much broader relaxation 
curve than that observed in the experiments. Our data 
indicate that a distribution of second order anisotropy 
with a non zero mode is needed to explain the MQT re- 
laxation in Mni2-ac. 

We want to note that we were able to fit the relaxation 
curves by using a discrete multi-peak distribution of E 
values similar to that expected from the solvent disor- 
der model However, it was necessary to include a 
Gaussian width to each peak of the distribution in order 
to fit of the data ^3 ■ The values of the peak centers of 
this distribution for resonance k = Q are Xd =-7.19 (- 
7.0525), xc,2 =-8.55 (-8.1749), and Xc,3 =-6.60 (-6.7995), 
(the values in parenthesis are extracted from ref. 16]). 
The width of each peak is Wi — Xc,i/50 and the height of 
each peak was taken to be proportional to the population 
of the corresponding isomer given in ref. ,16j . 

To conclude this subsection, we have shown that a 
multi-crossing LZ method allows the determination of 
the complete distribution of tunnel splittings for several 
resonances in SMMs. Our results suggest that a distri- 
bution of second order transverse anisotropics with a non 
zero mode is necessary to explain the experimental data. 
The solvent disorder model provides such a source and a 
discrete multi-peak distribution of tunnel splittings can 
be used to fit our data. However, an additional source 
of disorder (i.e. line-dislocations or point defects) that 
introduces a small broadening of the these peaks is also 
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necessary to model the experimental data. 



C. MQT Symmetry Measurements 

In this subsection we will present LZ relaxation exper- 
iments carried out in the pure quantum regime {T = 0.6 
K). In order to check the MQT symmetry imposed by the 
transverse terms of the Hamiltonian we have studied the 
LZ relaxation of the magnetization by sweeping an exter- 
nal longitudinal field, Hl, at a constant rate, a, across 
a resonance k in the presence of an external transverse 
field, Ht, applied at arbitrary directions, 0, with respect 
to the crystallographic axes of a Mni2-ac single crystal. 
As we have shown in section II (i.e. see fig. 4), MQT has 
an oscillatory response as a function of the orientation of 
a transverse field. This leads to maxima and minima in 
the MQT relaxation rates whose positions and symmetry 
depend on the transverse anisotropy term that generates 
the tunnel splittings. To recall, a fourth order anisotropy 
term would generate a fourfold rotation pattern in the 
MQT probability with maxima spaced by 90 degrees, 
while a twofold rotation pattern with spacing between 
maxima of 180 degrees is expected from a second or- 
der anisotropy term in the Hamiltonian. When two non- 
collinear different order transverse anisotropy terms are 
present, the symmetry of the MQT relaxation rates de- 
pends on the relative orientation between the anisotropy 
axes, as discussed in Sec. IIB. 

A single crystal of Mni2-ac SMMs was placed over a 
high sensitivity micro-Hall magnetometer as described 
in the first paragraph of this section. However, for these 
studies it is very important to know the exact orientation 
of the crystallographic axes with respect to the direction 
of the external magnetic field. Fig. 21 shows a sketch of 
the orientation of the c-axis of the crystal with respect 
to the external magnetic field. Note that, in Mni2-ac, 
the c-axis corresponds to the easy magnetic axis of the 
molecules. However, there is a misalignment of (p ^ 12° 
between the transverse magnetic axes (imposed by the 
fourth order anisotropy of the molecules) and the crys- 
tallographic axes. This is also shown in fig. 21. 

During the manual alignment of the crystal one of 
its faces was placed coplanar with the micro-Hall sen- 
sor plane with the help of a microscope. However, there 
exists an uncertainty of about ±5 degrees in this orien- 
tation. This misalignment is represented by the angles 
9 and cr in fig. 21 and can be determined experimen- 
tally through the magnetic measurements described be- 
low. The main implication of this misalignment is that 
there is a transverse field component due to the high lon- 
gitudinal field applied along the z-axis that depends on 
the angles 9, a and (p. The latter is the angle of ap- 
plication of the external transverse field with respect to 
the X-axis. Due to this, the total transverse field felt by 
the molecules in the presence of a constant transverse 
field, Ht, applied along and with a longitudinal field. 




FIG. 21: Schematic representation of the orientation of the 
c-axis (easy axis) of a crystal with respect to the axes of the 
applied external magnetic field {x,y,z). The misalignment 
is determined by 9 and a. The misalignment between the 
transverse magnetic axes of the Mni2-ac molecules, and one 
of the faces of the crystal {ip ~ 12°), is also shown. 

Hl = HzCos9 ^ Hz {cos9 ^ 1 tor 9 small) is given by, 
H^ - [HTCOs{(j) - cr) - HLsin9f + [HTsin{(t> - 

(14) 

Thus, the transverse field felt by the molecules will 
have a minimum value for (f> — cr and a maximum for 
(j) — cr-|- 180° with Hl > 0. The opposite situation would 
be found for Hl < 0. As we will show, once the mis- 
alignment angles are known, an algorithm can be used 
to correct the applied fields in order to have a constant 
transverse field during the measurements, independent of 
the angle of application of the external fields. 

In our first experiment we have studied the MQT re- 
laxation rates of several resonances in the presence of a 
constant transverse field applied at arbitrary directions 
with respect to the crystallographic axes. The exper- 
iment was done as follows: We start with the initial 
magnetization of the sample equal to positive satura- 
tion, Minitiai = +Ms, by applying a high longitudinal 
magnetic field, Hl — For this initial situation, all of 
the molecules of the crystal were in the m = -1-10 level 
in one of the energy wells. Then we turn on a transverse 
magnetic field, Ht = 0.4 T, appfied along a direction, 
0, and sweep the longitudinal field at a constant rate, 
a — 6.6 X 10^"^ T/s. We measured the magnetization 
change in several resonances k and determined the MQT 

probability, PlZ = {Mbefore - Mafter)/ {Mbefore " Meq), 

where, in this case, M^q = —Mg. Note here that all the 
molecules within the crystal contribute to the relaxation. 
We repeated this procedure for different angles from 
to 360 degrees. 

The behavior of the measured MQT probability as a 
function of (f) is shown in fig. 22 for resonances k = 
5 and fc = 6. The results clearly show fourfold max- 
ima in the tunneling probability spaced by 90 degrees 
(0,„„^ = 60°, 150°, 240° and 330°) for both resonances. 
Note that fc = 5 and 6 are the first observed resonances 
and ~ 35% of the magnetization relaxes. This means that 
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FIG. 22: Measured MQT relaxation probability for reso- 
nances k — 6 and fc = 7, as a function of the orientation 
of the applied transverse field, Ht = 0.4 T, relative to one of 
the faces of the crystal. The measurements where carried out 
starting from saturation (all the molecules within the crystal 
contributed to the relaxation) . 

molecules that contribute to this relaxation are mainly 
those with among the biggest tunnel splitting values 
within the distribution. 

There is also a one-fold contribution that is represented 
by a continuous line in fig. 22 with the result for reso- 
nance k — 6 44]. This oscillation is due to the mis- 
alignment of the c-axis of the crystal with respect to the 
applied magnetic field. This misalignment is represented 
by the angles 6 and cr and can be determined through 
magnetic measurements. From the results shown in fig. 
22 we find a = 80° (where we observe the maximum value 
of the one-fold contribution). To obtain 9, we measured 
the behavior of the MQT probability of resonance fc = 6 
for different magnitudes of the transverse field from -0.4 
to 0.4 T applied along cf) — a. The transverse field value 
for which the probability is minimum (null real trans- 
verse field) gives the angle of misalignment 6 through 
HT{Pmin) ~ HlsiuO, where HLsinO is the transverse 
projection of the longitudinal field. The value extracted 
for the misalignment is 6* = 0.3°. The effect of this mis- 
alignment on the transverse field is larger the greater the 
resonance number. We do not show the results for higher 
resonances, such as fc = 8, because the fourfold symmetry 
is almost unobservable due to the high one-fold contri- 
bution of the misalignment. 

In order to measure the response of the MQT prob- 
ability versus the angle (p in other parts of the tunnel 
splitting distribution we have conducted LZ relaxation 
experiments in the following manner. We select a small 
fraction of molecules with the smallest tunnel splittings 
of the distribution, in contrast to the biggest values that 
were analyzed in the previous experiments. The mea- 
surement method is presented in fig. 23. We start with 
M = —Ms by applying a high negative longitudinal mag- 
netic field, Hl = -5 T. Then we sweep the longitudi- 
nal field at a constant rate, a = 6.6 x 10~^ T/s, up to 
Hl = + 4.2 T (just after crossing resonance fc = 8) and 
sweep it back to zero, crossing again resonances fc = 8, 




FIG. 23: Selection of 30% of the molecules with the smallest 
tunnel splittings within the distribution. Starting at M = 
—Ms, the longitudinal field is swept from to 4.2 T, then back 
to in the absence of a transverse field (black line). After this 
procedure, only those molecules that have not relaxed (30%) 
remain in the metastable well. A transverse field of 0.4 T is 
then applied at an angle <j>, and the longitudinal field is swept 
again to a high positive value. The process is repeated for 
different (j) angles (lines with different colors). 



7, 6... The whole selection process is done in the absence 
of a transverse field. After that, the final magnetization 
of the sample is M = QAMg. This means that 70% of 
the molecules have relaxed to the stable well or, in other 
words, only 30% of the molecules have remained in the 
metastable well and will contribute to further relaxation. 
The latter are those molecules with the smallest tunnel 
splitting values within the distribution. After this, we 
turn on a transverse field, Ht = 0.4 T, applied at an 
arbitrary angle, 0, with respect to the crystallographic 
axes of the sample. Then we sweep the longitudinal field 
to a high positive value, crossing the resonances again. 
We repeated this procedure for different orientations of 
the transverse field from to 360 degrees. 

The MQT probability of resonance fc = 6 is shown as 
a function of in fig. 24. The results show the same 
fourfold symmetry pattern with maxima placed at the 
same positions as those in the experiment shown in fig. 
22. As these two experiments study the relaxation of 
two different parts of the distribution of tunnel splittings 
(low and high ends of the distribution) , we can conclude 
that the fourfold symmetry of MQT is a property of a 
significant fraction of the molecules within the crystal. 
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FIG. 24: MQT probability for resonance k = 7 versus the 
orientation of the applied transverse field, Ht = 0.4 T, rel- 
ative to one of the faces of the crystal. A previous selection 
process was used to study only 30% of the molecules with the 
smallest tunnel splitting values within the distribution. The 
results were extracted from longitudinal magnetic field relax- 
ation curves recorded at a constant sweep rate, a = 6.6 x 10~^ 
T/s. 

In principle, the four fourfold rotation pattern is con- 
sistent with a fourth order transverse anisotropy term, 
C{S^ -f 5*1), in the spin-Hamiltonian (see fig. 4b). For 
positive C, the four maxima generated by this term 
should be at (j)„,ax = 45°, 135°, 225° and 315°. There is 
a difference of 15 degrees between these values and those 
observed in our experiments, 4>max — 60°, 150°, 240° and 
330°. This is due to the misahgnment, (p — 12°, between 
the hard anisotropy axis of the molecules and one of the 
faces of the crystal, which we use as the origin of our 
rotation. There still is a difference of 3 degrees that is 
within the accuracy with which we orient the crystal (±5 
degrees). However, the value of C 3 x 10~^ K cannot 
explain the difference between the maximal and minimal 
magnitudes of the measured MQT probability. The ob- 
served normahzed changes, {Pmax - Pm'm)/Pmax, in the 
experiment with 100% of the molecules contributing to 
the relaxation are ~0.9 for k — 5 and ^0.6 for fc = 6 (see 
fig. 22). Whereas, with this value of C, we expect this 
change to be within the noise of the measurement. 

In order to determine whether this symmetry is in- 
trinsic to the Mni2-ac molecules and its origin, we have 
carried out experiments designed to select different parts 
of the distribution by using transverse fields in the se- 
lection process. Note that, in the previous experiment, 
we selected the SMMs with the smallest tunnel splittings 
in the absence of transverse fields. Now, we use a selec- 
tion transverse field (STF), Hstf, applied at an angle, 
(f'STF, during the preparation of the initial state of the 
system. In this case, those molecules with the medium 
anisotropy axis aligned with the STF have larger tunnel 
splitting values (larger relaxation probability) and can be 
selected for further study. 

The measurements with this selection process are 
shown in fig. 25. First we apply a high longitudinal field 
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FIG. 25: Selection processes carried out by sweeping a pos- 
itive longitudinal field across resonances fe = 5 and k — 6, 
in the presence of a selection transverse field, Hstf = 0.6 
T, applied at a selection angle 4>stf = 60°. Both 10% and 
50% of the molecules with the largest tunnel splitting values 
were selected in separate processes for this experiment. These 
populations were then studied on the negative side of the lon- 
gitudinal field hysteresis curve, in the presence of a transverse 
field, Ht ~ 0.3 T, applied at different angles (j> from to 360 
degrees. 



to saturate the magnetization of the system and sweep 
the field back to zero, having at the end a magnetization, 
Af = —Ms- Then we turn on a selection transverse field, 
Hstf = 0.6 T, applied at the angle (j)STF — 60°, where 
one of the four maxima were observed in the experiments 
carried out with the whole sample, and we sweep the lon- 
gitudinal field at a constant rate to a positive value and 
sweep back to zero. The value of this field is chosen de- 
pending on how much relaxation we want in the selection 
process: for a selection of 50% of the biggest splittings, 
we sweep the field up to 3.2 T allowing the system to re- 
lax in resonances fc = 5 and 6; for a selection of the 10% 
of the biggest splittings, we only allow the system to relax 
in resonance fc = 5. The final states of the magnetization 
are M = and M = -O.SMg respectively. After the 
selection process we sweep the longitudinal field down to 
-5.5 T at a constant rate, a = 6.6 x 10~^, in the presence 
of a transverse field, Ht = 0.3 T, applied at difference 
angles (j> with respect to one of the faces of the crystal. 
We repeated this procedure for different angles (p from 
to 360 degrees. In a separate experiment, we repeated 
the selection of 50% of the molecules with the largest 
splittings within the distribution by applying a selection 
transverse field, Hstf, along the angle in which a com- 
plementary maximum was observed, 4>stf = 150°. 

The results are shown in fig. 26. In all selections the 
MQT probability shows a twofold rotation pattern with 
maxima spaced by 180 degrees. For the selection in which 
Hstf is applied along (psTF = 60°, the two maxima are 
at (pmax,! — 60° and 4'max,2 — 240°, for both fractions 
of molecules selected (10% and 50%). Moreover, when 
the selection field is applied along the position of a com- 
plementary fourfold maxima, (psTF = 150°, the twofold 
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FIG. 26: Behavior of the MQT probabihty versus the orien- 
tation of a transverse field, Ht ~ 0.3 T. The molecules which 
relax were previously selected from the tunnel splitting distri- 
bution using a selection transverse field, Hstf, applied along 
the directions where the two complementary maxima where 
observed in the experiments on the whole sample (fig. 24), 
i.e. 4>sTF = 60° (upper figure, for both 10% and 50% of the 
largest splittings in the distribution), and 4>stf ~ 150° (lower 
figure, for 50% of the largest splittings). 



maxima, 4>max,i — 150° and 4'max,2 — 330°, are displaced 
by 90 degrees with respect to the previous case. 

The observation of a twofold rotation pattern in the 
MQT probability is clear evidence of a second order 
transverse anisotropy lower than that imposed by the 
site symmetry of the molecule (four- fold) . This is in ex- 
cellent agreement with the solvent disorder model pro- 
posed by Cornia et al. [l^ where there is an expectation 
of equal populations of molecules with opposite signs of 
the second order anisotropy. A considerable increase of 
the change of probability between maxima and minima 
is also observed as the fraction of molecules with highest 
tunnel splitting values becomes smaller, indicating the 
fact that the molecules with largest splitting values have 
bigger values of the second order anisotropy. 

In section II we stated that incommensurate anisotropy 
terms would considerably modify the magnetic response 
of the molecules to the orientation of a transverse mag- 
netic field, depending on the angle of misalignment, /3, 
between both anisotropics (see eqs. (5) and (6) and fig. 
5) . The symmetry of the MQT probability expected from 
this model depends on several parameters like the reso- 
nance number, fc, and/or the misalignment angle, /?. It 
turns out from this model that for \E\ < 30mK (bigger 
than the maximum E value expected from the solvent 
disorder model) the symmetry of the MQT probability 
mainly depends on the angle of misalignment, /3, go- 
ing from fourfold for small (3 values to twofold for big 
/3 values. Taking a given E value, the transition be- 



FIG. 27: Left: Experimentally determined angle dependence 
of the MQT probability for resonance fc = 6 for the two com- 
plementary directions of the transverse selection field. The 
one-fold contribution arising from the misalignment angles 9 
and a has been corrected for clarity. Right: Calculated tun- 
nel splitting for resonance k = 6, A^^s, for /? = —30° and 
opposite signs of E; the calculations assume the same values 
of Ht, D, E and C as those in fig. 9. 



tween fourfold to twofold maxima patterns as a function 
of (3 depends on the resonance number k. For exam- 
ple, for resonance fc = (fig. 8A), the fourfold max- 
ima pattern is observed for every misalignment value, 
where a slight twofold modulation is due to the second 
order anisotropy. For bigger resonances (i.e. A: = 6), 
the transition between these two MQT probability sym- 
metries is cleaner (see fig. 9). For misalignment angles, 
1/3 1 < 20°, the MQT probability shows fourfold symme- 
try modulated by the second order anisotropy. However, 
for angles \(3\ > 20°, the MQT probability symmetry is 
twofold. Consequently, the observation of twofold sym- 
metry in our experiments shows that the angle of mis- 
alignment between both anisotropy terms is greater than 
20 degrees. In section IV we will present high frequency 
EPR experiments that show that the angle of this mis- 
alignment is (3 = —30°. In fig. 27 we show a compar- 
ison between the experimental observation of resonance 
/c = 6 for both transverse selection fields with 50% of 
the biggest splittings (left polar plot) and the calculated 
splitting (right polar plot) corresponding to resonance 
k = 6 and /? — —30° with opposite signs of \E\ — 3mK. 
The difference of ~15 degrees between the experimental 
results and the calculations is also shown in this figure. 
As we said before, this difference is due to the fact that 
we measure the angle (p with respect to one of the faces 
of the crystal while the transverse magnetic axes of the 
molecules are rotated from the faces of the crystal by ^12 
degrees. 

An estimation of the values of E needed to explain the 
experimental observations of the oscillation of the MQT 
probability presented in this subsection are: a) E 0.5 
mK for the result corresponding to the 30% of the small- 
est splittings of the distribution (fig. 24), b) i? ~ 2.5 
mK for the result with the whole distribution (fig. 22, 
upper curve) and with 50% of the biggest splitting in 
the selection with Hstf = 60° and 50° (fig. 26), and 
c) ~ 10 mK for the result corresponding to 10% of 
the biggest splittings of the distribution (fig. 26, solid 
circles). These values are in excellent agreement with 
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FIG. 28: MQT probability for resonances A; = 5, 6 and 7 as a 
function of the magnitude of a transverse field applied along 
= 60° (open symbols) and <^ = 150° (solid symbols). The 
initial state of the sample was prepared by selection of 50% 
(left) and 10% (right) of the molecules with the largest tunnel 
splitting values, using a selection transverse field, Hstf ~ 
0.6 T, applied along (psTF = 60°. The two transverse field 
orientations correspond to the first maximum (0 = 60°) and 
minimum (0 — 150°) of the twofold rotation pattern. 

the results obtained by high frequency EPR experiments 
(presented in section IV) and by recent density functional 
theory calculations l35|. However, they are slightly larger 
than the values initially proposed in ref. [l6|. 



D. Berry Phase Measurements 

The ability to select a subset of molecules with a nar- 
row distribution of tunnel splittings and, for example, 
different signs of the second order transverse anisotropy 
allows us to study the behavior of the MQT relaxation as 
a function of a transverse field. As we anticipated in sec- 
tion II, quantum phase interference (Berry phase) would 
lead to zeros of the tunnel splittings (i.e. absence of mag- 
netic relaxation) for several values of a transverse field 
applied along the hard anisotropy axis of the molecules, 
and we discussed how the pattern of the oscillations can 
be modified by the presence of two non-coUinear trans- 
verse anisotropics. 

For these experiments we have used the same prepara- 
tions of the initial states of the system as those shown in 
fig. 25. These initial states correspond to a selection of 
50% and 10% of the molecules having the larger tunnel 
splitting values within the distribution. The selection of 
both initial states is done by applying a selection trans- 
verse field, Hstf — 0.6 T, at an angle 0stf = 60°. 
As we have shown, this procedure mainly selects those 
molecules with one i?-sign, as is observed in the twofold 
transverse field rotation pattern of the MQT probability 
of fig. 26, with maxima at (j)max = 60°, 240° and min- 



ima at (j)max = 150°, 330°. To study the behavior of the 
MQT probability as a function of the magnitude of the 
applied transverse field we apply a transverse field, Ht, 
along the direction of the first maximum, = 60°, after 
the selection process. We then sweep the longitudinal 
field at a constant rate to —5.5 T, measuring the change 
of magnetization in each resonance crossing. We follow 
the same procedure for different values of the transverse 
field from Ht = to 0.7 T. Moreover, we repeated the 
same measurement by applying the transverse field along 
the direction of the first minima tj) — 150° jS^. The re- 
sults of the MQT probability for resonances A: = 5, 6 and 
7 are shown in fig. 28. In the left graph of the figure 
are the results obtained with 50% of the molecules with 
the largest tunnel splitting values. In the right graph we 
show the results for 10% of the biggest splittings. The 
measurements with the transverse field applied along the 
first of the twofold maxima, (f) — 60°, are represented 
by open symbols, while the solid symbols correspond to 
measurements with the transverse field applied along the 
first of the minima, 6 = 150°. 



There are several important aspects to this figure, a) 
There is vertical shift between the curves correspond- 
ing to different directions of application of the transverse 
field. The shift is bigger in the case of the selection of 
the 10% of the biggest splittings within the distribution. 
This is consistent with the observations of fig. 26, and 
supports the assumption of a distribution in the mag- 
nitude of the second order anisotropy. b) MQT proba- 
bility increases exponentially with the magnitude of the 
transverse field. This is expected from the exponential 
dependence of the LZ probability on the tunnel splitting 
shown in eq. (8), and the power law dependence of the 
tunnel splittingon the magnitude of the transverse field 
(see eq. (3)) jsj- There are significant deviations from 
the exponential behavior in the right hand graphic for 
a transverse field applied along the first of the twofold 
minima, (j) = 150°. The largest deviations are observed 
at fields Hp{k = 5) ~ 0.45 T, Hp{k = 6) - 0.3 T and 
Hp{k — 7) 0.35 T. This is reminiscent of the Berry 
phase observed in Fes HI]- 



The results shown on the right graph of fig. 28 can 
be compared with the calculations of the tunnel splitting 
versus the transverse field shown in fig. 13 that were 
performed by taking into account two incommensurate 
anisotropy terms in the Hamiltonian (eqs. (5) and (6)). 
In fact, the values, £' = 10 — 15 niK, used in these cal- 
culations were chosen according to the results shown in 
fig. 28 and agree with the values extracted from high fre- 
quency EPR measurements (section IV). The agreement 
between theory and experiments is very good and consti- 
tutes the first evidence of quantum interference phenom- 
ena in a SMM system with incommensurate transverse 
anisotropics. 
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E. Summary of magnetic relaxation experiments 

We have shown in this section that LZ magnetic re- 
laxation experiments allow us to determine the complete 
distribution of tunnel splittings in Mni2-ac. The results 
obtained through a multi-crossing LZ method show that 
a distribution of second order anisotropics with non-zero 
mode is required in order to explain our data, such as 
in the solvent disorder model proposed by Cornia et al. 
|l6j |. LZ relaxation experiments carried out in the pres- 
ence of a transverse field applied at arbitrary directions 
with respect to the crystallographic axes of the sample 
enabled studies of the symmetry of the MQT probability. 
We have shown that the MQT probability has a general 
fourfold rotation pattern as a function of the orienta- 
tion of a transverse field. This is associated with equal 
populations of molecules with opposite signs of a second 
order transverse anisotropy. The LZ method allows the 
selection of a subset of molecules with different values 
of the tunnel splitting for further study. By applying a 
transverse field in the selection process, we can select a 
fraction of molecules in the sample with lower symmetry 
and with different signs of E. Using this selection proce- 
dure we have studied a small fraction of molecules with 
one sign of E and with the largest tunnel splitting values 
within the distribution. These show an unusual Berry 
phase phenomena for several transverse field values that 
does not lead to complete zeroes in the tunnel splitting. 
Our results on the symmetry of MQT can be explained in 
terms of incommensurate transverse anisotropics in the 
Hamiltonian that explain, among other things, why the 
observed Berry phase phenomena does not depend in any 
simple way on that expected from either anisotropy term 
alone. 



IV. HIGH FREQUENCY EPR EXPERIMENTS 

High frequency (40 to 200 GHz) single crystal Electron 
Paramagnetic Resonance (EPR) measurements were car- 
ried out using a millimeter-wave vector network analyzer 
(MVNA) and a high sensitivity cavity perturbation tech- 
nique; this instrumentation is described elsewhere |46| . 
Temperature control in the range from 2 K to room tem- 
perature was achieved using a variable-flow cryostat. The 
magnetic field was provided by a horizontal supercon- 
ducting split-pair magnet with vertical access, enabling 
angle dependent studies (<0.1° resolution) and approx- 
imate alignment of the single crystal ^2^. In order to 
make accurate comparisons with the magnetization stud- 
ies presented in the preceding sections, all data presented 
in this section were performed on a single deuterated 
crystal of Mni2-ac (d-Mni2-ac), having approximate di- 
mensions 1 X 0.3 X 0.3 mm^. This crystal was selected 
from a batch of samples which had previously been re- 
moved from the mother liquor and stored for 1 year in 
a refrigerator (at 5° Celsius) prior to the measurement. 
This particular batch was grown using standard methods 



[4l|, albeit in a completely independent synthesis from 
the samples used for the magnetization measurements. 
The sample was separately cooled under vacuum from 
room temperature at 5 K/minute in one of two orien- 
tations for field rotation in i) the x — y plane, and ii) 
a plane perpendicular to the x ~ y plane. In the for- 
mer case, the sample was mounted on the side wall of a 
cylindrical TEOll cavity (center frequency = 53.1 GHz, 
Q'--^20,000) with its easy axis parallel to the cavity axis 
such that the microwave field Hi was aligned parallel to 
the sample's easy axis (and, therefore, perpendicular to 
the applied DC field). In the latter case, the sample was 
mounted on the end plate of the same cavity, and DC 
field rotation was carried out for angles close to the x — y 
plane (within 15°), with the microwave Hi field again 
parallel to the sample's easy z-axis. Field sweeps were 
restricted to 6.6 T due to limitations of the split-pair 
magnet. As will become apparent, the data obtained for 
the d-Mni2-ac are in qualitative agreement with earlier 
published results obtained for the hydrogenated Mni2-ac 



A. Magnetic symmetry measurements in the 
high-field limit 

In the preceding sections, it has been shown how the 
Landau-Zener method may be applied to SMMs in or- 
der to determine very weak transverse terms in Eq. (1). 
Moreover, the Landau-Zener method allows one to select 
molecules, based on the tunneling rates of the different 
species. Subsequently, by performing angle dependent 
studies on each sub-species, one can deduce the under- 
lying symmetries of the dominant tunneling matrix ele- 
ments. While this method is extremely powerful, it is ev- 
ident from the discussion in section II.C that, for systems 
with multiple sources of transverse anisotropy (intrinsic 
and extrinsic), deconvolution of the different contribu- 
tions to the tunnel splittings can be problematic, i.e. in 
Figs. 6 to 9 it is seen that competing E and C terms 
results in a competition between the two-fold and four- 
fold symmetries of these interactions. The reason for this 
competition is that, at low-fields, the E and C terms 
operate in different high orders of perturbation theory, 
thereby resulting in a complicated interaction between 
the two perturbations. 

At first sight, it is not obvious how EPR experiments, 
conducted in the GHz frequency range, could shed new 
light on the nature of transverse interactions which man- 
ifest themselves as miniscule tunnel splittings of order 
10^ Hz at low-fields. However, the so-called "tunnel- 
splittings" are measured by the Landau-Zener method at 
low fields wherein the transverse terms operate in very 
high orders of perturbation theory within an rriz basis 
(see Eq. (3)), where the quantization axis is defined by 
the global easy-axis of the crystal. For high magnetic 
fields applied in the transverse direction, such a picture 
is no longer valid due to the conflicting symmetries im- 
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posed by the crystal field and the applied field. Herein 
lies the beauty of the high-field EPR technique. By ap- 
plying a sufficiently strong transverse field, one can reach 
a limit in which the appropriate basis of spin states is de- 
fined by a quantization axis parallel to the applied field, 
i.e. the x-direction. In such a limit, transverse zero-field 
interactions operate in zeroth-order. Consequently, their 
effects may be rather strong. As an illustration of this 
point, consider the simplest zero-field Hamiltonian: 

n = -DSl + E{Sl-S^)-gfiBH,Sx, (15) 

Making a substitution for in terms of {S^ + Sy), one 
obtains: 

n=^{D + 2,E)Sl + \{D ~ E){Sl ~ Si) - gfisH^S, , 

(16) 

which can be re-written as: 

n^^{D + 3E)Sl - g^iBH,S^ + H'r , (17) 

This equation is diagonal in Sx, and has the same form 
as Eq. (15). Therefore, to lowest order, the high-field 
eigenvalues will be given by: 

^{rrix) « ^(f -f ?>E)ml ~ guBH^m^ , (18) 

Similar arguments hold for higher order transverse terms. 
Thus, the transverse high-field EPR spectra provide per- 
haps the most direct means of measuring these transverse 
terms. Indeed, this represents one of the more illustra- 
tive examples of the importance of high-field EPR as a 
spectroscopic tool for studying quantum magnetism. Not 
only does the effect of the zero-field transverse terms shift 
to zeroth order, but the symmetry of such interactions is 
also preserved. This is best illustrated using the same 
example as above, with the field applied along the y-axis 
instead of the x-axis. In this case, Eq. (15) may be 
re- written: 

H=^{D-3E)Sl~g^iBHySy+n'i , (19) 

giving 

e('^y) ~ ~ iE)my ~ gfiBHyfUy , (20) 

Thus, from Eqs. (18) and (20), it is apparent that the 
influence of the (E) term changes sign upon rotating the 
applied field from the x-axis of the £'-tensor to the y- 
axis. This two-fold behavior is not unexpected; indeed, 
it is obtained also from the exact diagonalization cal- 
culations shown in Figs. 6 to 8, which correspond to 
the high-field/frequency limit discussed here. However, 
unlike lower-field calculations, the 2"'' and 4*'' order in- 
teractions decouple completely at high-fields, as do other 
transverse interactions. Thus, one may consider their ef- 
fects completely independently. This point is illustrated 



in Fig. 8, where it can be seen that the effect of an in- 
trinsic fourth-order anisotropy is to produce a ground-to- 
first-excited-state splitting (Afe=o) which oscillates as a 
function of the field orientation within the hard/medium 
plane, with a periodicity of 90°. The disorder- induced 
rhombic anisotropy, meanwhile, has no affect on this four- 
fold behavior. It simply causes a two- fold modulation of 
the tunnel splitting, which superimposes onto the four- 
fold behavior caused by the fourth-order interaction. Be- 
cause of the complete independence of these effects, one 
can in principle determine any misalignment, /?, between 
the hard axes associated with the second and fourth order 
transverse anisotropies, as illustrated in Fig. 8. 

In order to make direct com pari sons with other spec- 
troscopic studies (e.g. neutron [l^ and EPR '7,15,1^) 
we re- write the Hamiltonian of Eq. (1) in the following 
form: 

n = D'[Sl - ^SiS + 1)] + EOl' + Biol + B|6| + Hz 

(21) 

where Tiz = ^jlbB ■ g ■ S is the Zeeman term due to 
the applied magnetic field, EO2 is the first term in 
eq. (5), and are the Steven's operators, of order B 
in the spin operators and possess A-fold symmetry (i.e. 
O4 = \{S\_ -t- S\)). The uniaxial parameter D' is not 
exactly the same as the D parameter in Eq. (1) [or Eqs. 
(15)-(20)]. This is due to the occurrence of an Si term in 
the 0\ Steven's operator, i.e. the presence of a significant 
fourth-order axial anisotropy has the effect of renormal- 
izing the quadratic m-dependence of the barrier [iR . 

Fig. 29 displays the raw microwave absorption ob- 
tained for different field orientations within the hard 
{x—y) plane of the sample; the temperature was 15 K and 
the frequency was 51.3 GHz in every case. The peaks in 
absorption correspond to EPR. The data were obtained 
at 7.5° intervals, where the angle (f) refers to the field 
orientation relative to one of the flat edges of the square 
cross section of the sample. The resonances have been la- 
beled according to the scheme described in ref. [23| . For 
fields applied approximately parallel to the hard plane, 
only a— resonances are observed (/?— resonances appear 
for field rotation away from the hard plane, see below). 
The highest field peak, Q!8, corresponds to an excitation 
between levels which evolve from the zero field = ±9 
zero-field doublet. The transition from the ground state, 
rriz = ±10 (fc = resonance), is not observed within the 
available field range for these experiments; at 51 GHz, 
its expected position is at Ht ~9 T. For a detailed dis- 
cussion of the resonance labeling scheme, as well as the 
temperature, frequency, field and field orientation depen- 
dence of the EPR spectra for Mni2-ac, refer to |27j . 

Immediately apparent from fig. 29 is a four-fold 
variation in the positions of each cluster of resonances 
(a8, a6, etc.). Note that each of the resonances ex- 
hibit fine structures, which also depend on the field 
orientation 0; these are related to the disorder in the 
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FIG. 29: Microwave absorption obtained for different field 
orientations within the hard (xy-) plane of the sample; the 
temperature was 15 K and the frequency was 51.3 GHz in 
every case. The peaks in absorption correspond to EPR. The 
data were obtained at 7.5° intervals, where the angle (p refers 
to the field orientation relative to one of the fiat edges of the 
square cross section of the sample. The resonances have been 
labeled according to the scheme described in ref. The red 
traces correspond to field orientations approximately parallel 
to the hard/medium axes of the E tensor, i.e. orientations 
corresponding to the maximum splitting of the low and high 
field shoulders. The blue traces correspond to orientations of 
the hard axes of the B2 tensor. 



crystal, which we discuss further below. The four-fold 
shifts are due to the intrinsic fourth-order transverse 
anisotropy B4O4, as has previously been established for 
h-Mni2-ac ,19]. Fig. 30 shows a color contour plot of the 
absorption intensity versus magnetic field strength and 
the azimuthal angle (p. Superimposed on the absorption 
maxima are two kinds of fit to the ^-dependence of 
each peak. The solid blue curves were obtained simply 
by fitting the positions of the central peaks in Fig. 29 
with pure sine functions having four-fold periodicity. 
The fits represented by horizontal bars, meanwhile, were 
obtained by exact diagonalization of the Hamiltonian 
matrix, assuming accepted values for the zero-field 
parameters D' and {D' = -0.455 cm~^ and 5° = 
2xl0~^ cm~^). The fourth-order transverse parameter 
B4 = 3.2 X 10~^ cm~^ is the only unknown parameter in 
the fit. D' and were verified independently from easy 
axis measurements, and all peak positions are consistent 
with a single value of Bf. We note that this value is in 
precise agreement with that found for h-Mni2-ac, as is 
to be expected for this fourth-order interaction which 
is related to the intrinsic symmetry of the Mni20i2 
molecule. From the maxima and minima in the peak 
position shifts induced by the BfOl term, we estimate 
that the hard and medium directions of this intrinsic 
crystal field interaction are oriented at 4>hc = —4.5° ±5° 
{+i90°) and 0AfC = 41.5° ± 5° {+i90°) relative to the 
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FIG. 30: Color contour plot of the absorption intensity (see 
Fig. 29) versus magnetic field strength and the azimuthal 
angle 0; the darker shades correspond to stronger absorption. 
Superimposed on the absorption maxima are fits to the 0- 
dependence of the central positions of each peak (solid blue 
lines), as well as fits to the positions of the shoulders (solid 
red lines). The horizontal bars are fits to Eq. (21) (see main 
text for explanation). The approximate orientations of the 
hard axes corresponding to the E (HE) and Bf (HC) tensors 
are indicated. The open circles represent recent data points 
obtained for h-Mni2-ac li^. 



square edges of a typical single crystal sample. 

Next we turn to the angle-dependent fine structures 
which are very apparent in the ranges — 300° — 330° 
and (/) = 30° - 60° in fig. 29. The first point to note is 
the fact that we see shoulders on both the high and low- 
field sides of the main peaks in these angle ranges. This 
contrasts the situation found from our earlier studies of 
h-Mni2-ac I!) . where only high-field shoulders were ob- 
served. We comment on these differences at the end of 
this section. We first discuss the origin of the angle de- 
pendence of the shoulders, which are very apparent in 
the ranges (j) = 300° - 330° and = 30° - 60° in fig. 29; 
a more in-dcpth discussion can be found in refs. [19|, |27| . 
The high- and low-field shoulders are due to the n ~ 1 
and n = 3 Mni2-ac hydrogen-bonding variants in Cor- 
nia's solvent disorder model [l^. These variants, which 
comprise roughly 50% of the total molecules, are thought 
to possess a significant rhombic anisotropy due to the 
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reduced symmetry of the surrounding hydrogen bonded 
acetic acid solvent molecules. The second-order operator, 
O2, associated with the rhombic distortion gives rise to 
two-fold behavior, as was clearly demonstrated from the 
magnetic measurements in the previous section of this 
paper. It is important to recognize that the disorder is 
discrete, since the acetic acid can only bond at four po- 
sitions on the Mni2 molecule. Thus, one expects only 
two hard directions associated with the n = 1 and n = 3 
low symmetry variants, which are separated by 90°. In 
HFEPR, the 0| operator causes shifts in the peak po- 
sitions. When the applied transverse field is parallel to 
one of the hard axes, it is obviously perpendicular to 
the other, which causes shifts in the HFEPR intensity to 
both the low- and high- field sides of the central peaks, 
hence the shoulders. When the field is applied in be- 
tween these two directions, the EPR intensity due to the 
low-symmetry variants collapses into the central portion 
of the peak, hence the disappearance of the shoulders 
every 90° (see figs. 29 and 30). In actual fact, this ap- 
parent four-fold behavior refiects the two-fold nature of 
the rhombic distortion caused by the hydrogen bonding 
acetic acid molecules, with the EPR intensity for a given 
variant shifting from the low (high) to the high (low) field 
side of the main peak every 90°, i.e. the periodicity is ac- 
tually 180°. This two-fold behavior is then superimposed 
on the intrinsic four-fold periodicity, as can be seen from 
the red curves in fig. 30. 

An important point to note from the hard-plane rota- 
tion data is the phase shift between the four-fold modu- 
lation of the central peak position (blue curves in fig. 30) 
and the two-fold shifts of the low and high field shoul- 
ders (red curves in fig. 30). Indeed, this is one of the 
main points of this article. As discussed in section II. C, 
this difference indicates a ~ 27° ± 3° mis-alignment of 
the 0| and Of tensors, as was originally suspected for 
h-Mni2-ac 0, . The present study provides further 
support for this finding, thereby illustrating the remark- 
able differences between the global and local symmetries 
of Mni2-ac. 

In order to make quantitative comparisons between 
the disorder-induced effects in h-Mni2-ac and d-Mni2- 
ac, we performed a single fit to each of the peaks in 
fig. 30 via exact diagonalization of the Hamiltonian ma- 
trix [Eq. (21)]. This fit is represented by the horizon- 
tal bars in fig. 30. Our procedure obviously takes into 
account the mis- alignments of the O2 and O4 tensors, 
as described in section II-C. Thus, the employed Hamil- 
tonian is subtly different from the standard form used 
by most spectroscopists, which may explain slight dif- 
ferences in the obtained Hamiltonian parameters. The 
only free parameters in the fit were then the E and Bf 
coefficients (corresponding to the O2 and Of tensors), 
for which we obtained the values ±0.014(2) cm~^ and 
±3.2(5) X 10~^ cm~^ respectively. The E value is sig- 
nificantly l arg er than the one obtained in earlier experi- 
ments ref. ri9]). We attribute some of this difference to 
the modified Hamiltonian used in our more recent fits, 
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FIG. 31: a) Hard-plane angle (</!>) dependence of the splitting 
of the high and low-field shoulders, for the a8 peak. This fig- 
ure may be compared with Fig. 2 of ref. |13|. which displays 
similar data for the h-Mni2-ac complex. In (a), the positions 
of the low- (black) and high-field (green) shoulders on a8 
have been separately fit with sine functions. The difference 
between these fits is displayed in (c), together with similar 
curves generated by the same procedure for peaks q6 and q4. 
The splittings plotted in (c) are a measure of the shifts caused 
by the disorder-induced rhombic term for the low-symmetry 
Cornia variants 

(see also Fig. 8); all curves are m agree- 
ment (±2°) as to the orientation of the hard/medium two- 
fold directions. In (b), the data in (a) are plotted in such a 
fashion as to illustrate the real two-fold nature of the angle 
dependence of the high- and low-field shoulders. The approx- 
imate orientation (0 — —4.5°) of one of the hard four-fold 
axes (HC) is indicated in (b), and the approximate orienta- 
tion {(j) = —31.5°) of one of the hard/medium two-fold axes 
(HE) is indicated in (c). 



which takes into account the mis- alignments of the O2 
and 0| tensors. However, much of the difference appears 
to be real. The larger .E-value found from the present 
study is somewhat surprising. However, as will be seen 
below, more recent measurements on a very fresh h-Mni2- 
ac sample are in excellent agreement with the value of 
±0.014(2) cm"-'^ found from this study (open blue circles 
in fig. 30). Thus, the difference is likely related to sam- 
ple quality and/or solvent loss. A tell-tale sign of the 
higher sample quality is the observation of both high- 
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and low-field shoulders on the main EPR peaks. In con- 
trast, only high-field shoulders were seen in the earlier 
experiments on h- Mni2-ac ^^H^- In fact, the absence 
of a low-field shoulder is discussed at some length in ref. 
|27|. where it is shown that this peak is unresolved from 
the broad low-field tail associated with the central por- 
tion of the peak. The reason for the asymmetry between 
the high and low-field shoulder is related to easy axis 
tilting caused by the solvent disorder (discussed in the 
next sub-section). Although both shoulders are seen in 
the present study, a clear asymmetry can be seen from 
the data in fig. 29. We suspect that the samples used 
in the earlier experiments may have suffered significant 
solvent loss, either upon cooling from room temperature 
under high vacuum, or simply as a result of being stored 
in air for more than 1 year prior to the measurements. 
Indeed, variations in the Z?-strain measured in different 
Mni2-ac samples has previously been reported by us [2ll |. 
Increased D-strain leads to broader EPR lines, thus prob- 
ably explaining why the shoulders are clearly resolved in 
the present investigation, but not in refs. [13; 123 ■ With- 
out two well resolved shoulders, it is likely that the E- 
strain was under- estimated in ref. or it could sim- 
ply be that the i?-strain is weaker in samples that suffer 
significant solvent loss. Based on more recent experi- 
ence working with SMMs containing considerably more 
volatile solvents, we have recently developed sample han- 
dling procedures which minimize solvent loss, e.g. encap- 
sulating samples in oil prior to cooling under atmospheric 
helium gas. The differences between the present measure- 
ments and those reported in refs. 0, |23| highlight the 
importance of sample handling. Indeed, it is likely that 
Mni2-ac samples prepared by different groups, and stud- 
ied by different techniques, exhibit significant differences 
in their solvent content, resulting in subtly different con- 
clusions concerning the quantum dynamics. Solvent loss 
probably also provides an explanation for the widths of 
the distributions of tunnel splittings found in the mag- 
netic relaxation experiments described in section III. For 
this reason, it will be advantageous to prepare Mni2-ac 
SMMs which do not exhibit such a dramatic dependence 
on solvent content. 



B. Easy-axis tilting 

In this section, we briefly present the results of mea- 
surements for fleld rotations away from the hard plane 
in order to illustrate the presence of a small distribu- 
tion of tilts of the easy axes of the d-Mni2-ac molecules. 
In the following figures, the polar angle 9 represents the 
angle between the applied DC magnetic field and the 
global easy axis of the crystal, i.e. 6 — 90° indicates the 
hard/medium-plane direction. Field rotation was per- 
formed in a plane approximately parallel to one of the 
large fiat surfaces of the needle-shaped crystal (i.e. = 0, 
or 90°, etc.). A more extensive discussion of the analysis 
and interpretation of the results of similar experiments 
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FIG. 32: Microwave absorption spectra obtained at 0.2° in- 
tervals over the range from 9 = 90° to 97° {(p ~ 5°); the 
temperature was 15 K and the frequency was 51.3 GHz in 
every case, and the traces are offset for clarity. The peaks in 
absorption correspond to EPR, and the labeling is discussed 
in ref. |2^. See main text for discussion of the data. 



for the h-Mni2-ac complex have been published elsewhere 
[27l |. Consequently, we present only data for d-Mni2- 
ac in this article, and discuss the implications without 
including detailed simulations, which will be published 
elsewhere |47l |. 

Fig. 32 shows a series of absorption spectra obtained at 
0.2° intervals over the range from 9 = 90° to 97°. Again, 
the peaks in absorption correspond to EPR, and the la- 
befing is discussed in ref. |27|. The quality of the data 
is noticeably poorer than the data in fig. 29, and is due 
to the positioning of the sample at a different location 
in the cavity, where the geometry of the electromagnetic 
fields are not quite optimal for EPR. Fig. 33 displays (a) 
a gray scale contour map and (b) a 3D plot represent- 
ing the same data shown in fig. 32 (the darker colored 
regions correspond to stronger EPR absorption), albeit 
for a wider range of angles (up to 9 = 105°). As dis- 
cussed at great length in ref. |23|, the transverse-field 
EPR spectra for Mni2-ac exhibit unusual selection rules 
for frequencies below 90 GHz: two series of resonances 
are separately observed (labeled a and (3) as the magnetic 
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FIG. 33: (a) A color contour map and (b) a 3D color surface 
plot representing the same data shown in Fig. 32 (the darker 
colored regions correspond to stronger EPR absorption). We 
plot the data in this manner for direct comparison with simu- 
lations shown in Fig. 6. The significant overlap of the a and 
l3 resonances suggests tilting of the molecules. See main text 
and ref. [2711 for discussion. 



field is tilted away from the hard plane. These selection 
rules are extremely sensitive to the field orientation for 
angles close to the hard plane. Thus, deviations from 
the behavior predicted by the giant spin model (Eq. (1)) 
provide evidence for tilts of the molecules. 

51.3 GHz simulations of the EPR spectra for an ide- 
alized Mni2-ac sample are displayed in fig. 34 — both 
a color contour map for the full 9 ~ 90° to 105° range 
(fig. 34.a), and a 3D view for the 9 = 90° to 97° range 
(fig. 34. b). These 15 K simulations were generated us- 
ing the accepted Hamiltonian parameters given above. 
A Gaussian line shape was employed with a line width 
typical for a single solvent-disorder variant in Cornia's 
model (see ref. (13 )■ While this line shape/width is 
clearly quite different from the experimental data, which 



is complicated by the contributions of several Mni2-ac 
variants, we note that the choice of line width/shape 
does not in any way affect the following analysis. The 
simulations indicate a range of about 1.6° (from 91.8° 
to 93.4°) over which the EPR intensity associated with 
both the a8 and a7 transitions is negligible. Thus, for a 
perfect crystal, without any disorder, one should expect 
a similar behavior in the actual EPR spectra. However, 
careful inspection of fig. 32 and 33 indicates a signifi- 
cant overlap of the a8 and /37 peaks in the 92° to 93° 
range. These two facts point to a spread in the orien- 
tations of the magnetic axes of the molecules, with a 
cut-off of at ~ 1.3° away from the global directions, i.e. 
we predict that, on average, the magnetic easy axes of 
the low-symmetry (disordered) molecules are tilted away 
from the crystallographic z-axis, and that the distribu- 
tion extends roughly 1.3°. Once again, a very similar 
behavior has been observed for h-Mni2-ac in ref. (27l| . 
albeit that the distribution extends to ~ 1.7°. We note 
that the h-Mui 9-ac experiments were conducted using a 
rotating cavity |50| to higher magnetic fields, where the 
experimental deviation from the simulation is even more 
pronounced. Thus, easy axis tilting appears to be a gen- 
eral feature in Mni2-ac and, as discussed earlier in this 
article, this likely provides an explanation for the lack of 
selection rules in the magnetization steps observed from 
hysteresis experiments. The slightly broader tilt distri- 
bution for h-Mni2-ac may be related to the increased 
solvent loss. 

Finally, we comment on the correlation between the 
disorder-induced rhombic anisotropy and the easy axis 
tilting. It has been shown by Park et al. [s^, that 
the magnetic anisotropy tensor absolutely determines the 
orientations of the principal magnetic axes. In other 
words, one expects the effects of the solvent disorder to 
be accompanied by local easy-axis tilting. This was also 
pointed out by Cornia et al. in their X-ray analysis 
of Mni2-ac. We have subsequently shown this to be the 
case from EPR studies of h-Mni2-ac 27], where it was 
shown that the high-field fine structures in the transverse 
EPR spectra displayed a different angle dependence com- 
pared to the main peaks. These findings suggested that 
each solvent disorder variant has a distinct angle depen- 
dence and, therefore, a distinct tilting behavior, i.e. the 
tilting and the anisotropy are correlated. This is again 
apparent for d-Mni2-ac from fig. 33. a, where the shoul- 
ders are visible as narrow horizontal streaks on the high 
field sides of each of the main peaks. It is quite evident 
that the narrow streaks (i.e. the shoulders) span a nar- 
rower angle range compared to the main peaks. In fact, 
the high field shoulders on the ag and 017 peaks exhibit 
a considerable range where neither is observed. Indeed, 
this range corresponds almost exactly to the 1.6° found 
from the simulations (fig. 34). It is important to note 
that this local tilting of the magnetic axes (caused by 
local symmetry lowering), is quite different from physi- 
cal tilting of the molecules caused e.g. by strains in the 
sample. 
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FIG. 34: (a) A color contour map and (b) a 3D color surface 
plot representing simulations of the data in Fig. 33, assuming 
that all molecules are aligned, i.e. no tilting. The simulations 
were generated using the Hamiltonian parameters given in the 
text, and the temperature and frequency are 15 K and 51.3 
GHz respectively. 

Comparisons between Figs. 32, 33 and 34 indicate 
that the angle dependence of the high-field shoulders 
agrees very well with the simulations. The reason 
for these differences between the main peaks and the 
shoulders is explained in ref. as being due to discrete 
easy axis tilting, wherein the tilting is confined to two 
orthogonal planes defined by the hard and medium 
directions of the associated disorder-induced Ol{E) 
zero-field tensor. Essentially, the high-field shoulders on 
the EPR peaks are due to molecules which are tilted in a 



plane which is approximately perpendicular to the plane 
of rotation of the applied magnetic field. Consequently, 
the tilts due to these molecules do not project onto the 
field rotation plane, i.e. the experiment is insensitive 
to tilts in this direction. Meanwhile, the low-field tail 
of the EPR peaks is due to molecules which tilt in the 
orthogonal plane. For these molecules, the tilts have a 
maximum projection on to the field rotation plane, i.e. 
the experiment is maximally sensitive to tilts in this 
direction. The present studies support the findings of 
the original study for h-Mni2-ac |22|, providing further 
confirmation for the discrete tilting idea. 



V. CONCLUSIONS 

In summary, these experiments provide a compre- 
hensive understanding of the factors that influence 
the symmetry of MQT in Mni2-ac, the first and most 
widely studied SMM. Interestingly, the reason so much 
attention has focused on Mni2 has been the high global 
symmetry of single crystals, as most known SMMs have 
lower site symmetry. The data presented here show 
that disorder lowers the symmetry locally and leads to 
an intricate interaction between transverse anisotropy 
terms: the first associated with disorder and the second 
intrinsic to an "ideal" Mni2-molecule. While at first 
sight a nuisance, these complexities have been interesting 
from many perspectives. From experiment, they have 
been a challenge to understand and characterize, and 
have required combining unique advanced and sensitive 
magnetic techniques that have been developed by the 
authors over many years. In particular, this research 
involved combining single crystal high-frequency and 
high-field EPR and low-temperature magnetometry, 
both with arbitrarily directed applied magnetic fields. 
The results obtained by magnetic measurements of MQT 
have had implications for EPR studies and vice-versa. 
From a theoretical perspective, given a structural 
model of the molecule and solvent environment, density 
functional theory has been able to capture many of 
the features observed in these experiments, including 
the magnitudes and form of the transverse magnetic 
anisotropics and the easy axis tilts [35l |. This includes 
the angle between the 2nd order and 4th order transverse 
anisotropics, which has been central to understanding 
the combined data set, as well as the easy axis tilts. 
This research thus represents an important milestone in 
our understanding of the factors that infiuence MQT in 
SMMs. 
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